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PREFACE. 


In physics and technics one is often faced with the problem 
of studying the periodical steady state of a mechanical or electrical 
system, corresponding to a sinusoidal disturbance attacking it. 
If the system is governed by linear differential equations with 
coefficients, not depending upon time, one can, as we kiiqw, 
by introducing the representative complex vectors of the time 
functions obtain new, simpler equations, which do not contain 
the variable of time. This mode of proceeding has, especially in^ 
electrotechnics, got an extraordinary practical importance, beingi 
the basis of thp^ so-called symbolic method in the alternating- 
current theory. 

It seems to be 0. Heaviside who first systematically forked 
out methods of treating non- periodical and transient phenoH^ena 
in a symbolic way. In later years Heaviside’s so-called operational 
calculus has beei\ subject to an increasing " interest, and several 
^tempts have been made to found the many operations occurring 
upon a more general and mathematically more satisfactory basis. 
Hence, the literature about operational calculus is nowadays rather 
copious. A selected collecfion of textbooks and recent treatises 
on the subject is given at the end of this book. 

The aim of the present work is to give an account of certain 
methods of integration, nearly related to the Heaviside operational 
calculus but more consequently building upon the symbolic 
method of the alternating-current theory. The author has tried 
to give a relatively simple but mathematically well-founded 
generalization of the principal ideas and notions of this method. 
Of course, the Fourier integral theorem offers a natural way to 
gain this aim, but, like D. P. Dalzell, among others, the author 
has found the modification of this theorem due to Mellin more 
adequate for the purpose, and thus based the development of the 
methods upon that form of the theorem. This will cause ordi- 
nary analytical methods of the theory of functions of a complex. 
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variable to play a prominent part in the treatment. This will 
also be the case with some procedures of the functional calculus, 
because there exists a near rel^itionship between operational 
methods and this branch of analysis. Indeed, as outlined by P. 
L4vy, it is possible to build up most of Heaviside’s operational 
calculus upon the functional calculus as a basis. 

In the course of the work the author has been influenced by 
many other treatises, to mention only the wellknown textbooks 
by Jeffreys, Carson, Bush and Berg. As to the practical applications 
these books and the present one complete one another in many 
respects. It ought to be noticed that in the following treatise the 
main interest concerns the methodical and mathematical side of 
the subject; therefore a more detailed discussion of the solutions 
and especially of their physical contents is generally let aside. 
It has been my intention to present many results under such 
analytic forms as usually occur in other textbooks, and the reader 
is once for all referred to standard works such as those by Riemann- 
Weber, Webster and Courant- Hilbert besides the books already 
mentioned. 

When trying to develop consequently the operational calculus 
as an extension of the symbolic method the author got reason 
to introduce the notion of integrator, essentiaPy corresponding to 
the representative complex vector of a sinusoidal function; this 
brings also with it a certain difference from other textbooks in 
the mode of writing operational equations. However, also the 
notion of operator has been frequently used, mostly signifying a 
linear functional operation, and therefore it seems justified to 
employ the name integrational and operational calculus on the 
methods to be presented. 

Some proofs and supplements of a more completing or special 
nature have been transferred to an Appendix at the end of the text. 

This course is essentially meant for electrical engineers that do 
not content themselves with a purely formal management of the 
operational calculus. Of course, it does not claim at all upon 
completeness, but if, when studied at the .side of other similar 
books, it could give the reader some new and useful aspects on 
the subject, its aim would be gained sufficiently well. 

Orebro, Aug. 1935. 


K. D. 



A C9urse of 

Integrational and Operational Calculus 
* with applications to problems of 
Physics and Electrotechnics. 


1. The integral theorem. General points of view. 

In the following pages we shall have to deal with the im- 
portant problem of calculating the effects of arbitrary external 
forces attacking given mechanical or electrical systems. We 
assume the dynamical processes to be governed by differential 
equations that are linear and have constant coefficients with 
respect to the variable of time t. The physical sense of the 
linearity is the 'validity of the principle of superposition, 
expressing that the effects of different forces will independently 
be added to one another. Hence our problem essentially consists 
in computing tl|,e effect of a single force attacking the system 
at any moment; the influence of certain prescribed initial con- 
ditions may also be included in this problem. Nothing prevents 
us from taking still one step towards a simplification, assuming 
the external force to keep'^lhe constant value 1; indeed, forces of 
any nature may easily be built up by elements of that type, 
represented by tbe discontinuous unit function or »iump func- 
tion)) H{t). 

Accordingly, we are going to work with a general class of 
functions possessing the characteristic property of being zero iden- 
tically up to a certain moment, say <— 0, when they change, in 
general discontinuously, into functions more or less regular for 
positive values of t. 

The practical methods of integration depend upon the possi- 
bility of algebraizing the expressions in the time-derivatives and 
-integrals. In that way, instead of ordinary differential equa- 
tions in t we get purely algebraical equations, and partial dif- 
ferential equations in t and other variables x, y, z are re- 
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placed by simpler equations not containing t. To permit the 
algebraization the time function in question must be decomposed 
into a number of exponential elements. Generally the external 
forces will be of such a nature that the classical Fourier inte- 
gral theorem can not be applied, but a variant af this theorem, 
due to Mellin, will prove to be very suitable for the purpose. 

The Fourier-Mellin theorem can be formulated in the following 
way, adjusted with respect to our special state of things (Cf. 
Appendix 1): 

Proposed the real function / {t) possesses the following pro- 
perties: 

1. f {t) and its first derivative are piecewise continuous in 
each finite interval. 

2. f (t) = i [f {t+0)+f {t — 0)] at a point of discontinuity. 


+ 00 

/ l/WI- 


e-'^^ dc exists, when oci < a < 


Then / {t) may be represented by the complex integral 


where 


a—i oo 


4- oo 


F{p) = p- / f (r) er-J>' dc. 


(1) 

( 2 ) 


In (1) we assume Ihe path of integration to be a straight line 
parallel to the imaginary axis of the p-plane. This line, how- 
ever, can be arbitrarily removed within the domain 

«x < Be (p) < as 

without changing the value of the integral; it may also be re- 
placed by any curve in this band between the extremities 
a ± t oo (Cf. fig. 1 ). 

Without essentially limiting the practical applicability of our 
methods we may generally presume more special properties of the 
function f(f) than those required by the integral theorem. Hence, 
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Fig. 1. 

as already mentioned, we can assume that f(f)=zO, when / < 0, 
and accordingly put the lower limit in (2) equal to zero. Further, 
for positive values of t we will prescribe f(t) to be continuous 
and to have a generally continuous and exponentially limited 
derivative or sometimes such derivatives of still higher orders. 

A function (p (/) is said to be exponentially limited when / + oo, 

if the inequality | (/) I < holds for all t > t .^, the real ex- 

ponent a being independent of t. 

To simplify our mode of expression we say that a function 
/(/) belongs to class n if, for positive values of /, it is continuous 
and possesses continuous derivatives up to the order (n — 1) and 
an at least piecewise continuous and exponentially limited n:th 
derivative. Then, as is readily proved, there must exist a num- 
ber a and a number /:!« such that 

I (0 I < when />/:}; (3) 

(* = 0, 1, . . n) 

With the notation /(O) we always understand the limit value 
of f(e) when + 0. 
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A function 4s said to belong to the class of order infiniiij if it 
has an unlimited sequence of continuous, uniformly bounded, and 
uniformly exponentially limited derivatives. 

The latter type of functions is most common in the appli- 
cations. When not otherwise stated we will presume the time- 
functions appearing in the following to belong to class nr 1 
at least. 

On account of this statement the interval ai < o < a 2 in the 
integral theorem may be replaced by an infinite interval :o > a. 
Here a has the same meaning as in (3). Hence, the line of in- 
tegration L (Cf. fig. 1) can be situated anywhere to the right of 
the line Re {p} == a. 

It is easily verified (Cf. Appendix 3) that a function F(p), 
being defined in accordance with (2) by any function /(O, has 
the following general properties: 

1. Within the half-plane Re {p) > a the function — — is 

P 

holomorphical. 

2. Within the half-plane Re {p) > Oo > a the function — 
tends uniformly to zero, whenlpj-^oo. 

Some simple examples may illustrate the» subdivision of a 
function / (i) into exponential elements by the integral theorem. 


A. The unit-junction H (t). 

This fundamental function is the simplest in the class of func- 
tions in question. It is defined by: 


H(t) = 


0 

a 


1 


<< 0 . 

t=0.) 


t>0. 


The value for t=0, expressing a regular discontinuity at the 
origin, is immaterial for us. According to our agreement above 
we always put ff (0)=1. 

For this function we get by (2): 


C 

F{p)=pJ' i 
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if <J > 0, and accordingly, by (1): 

G+i OO 

<*> 

G—i OO 

for any line of integration on the right of the imaginary axis. 

If (p (t) denotes a function defined for any real value of t, 
e. g. by an analytic expression, we get, when multiplying by the 
unit-function, a new function, (p{t)-H(t), which is zero for all 
negative values of t and coincides with q) (t) for positive values 
of L 


B. The exponential function e^^ H {t) (a real). 
One gets the formulas: 


F 


QC 

(p)=p.y^ . 


^(a— -P)* 


p — a 


G+l OO 


2m J p — ap 


(5) 


provided that a>a. 

Note. When « is negative, we also have a=0>«; hence the imaginary axis 
may be used as a path of integration, which means that the classical Fourier 
theorem is applicable. 


C. The functions sintot-H (t) and cosot-H (t). 

If o>0 we get from the integral theorem 

O’+i OO 

(6) 

G—i OO 


cos bit-H (t) 


(T+t OO 
>2 


^J_. /l^ 

2m J p*H- 




G—i OO 


H-w® p 


■dp. 


(7) 


It is easily verified in the different cases that the complex 
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integral represents the explicit time function on the left hand 
side. 

Evaluating complex integrals of the t 3 ^pe (1) is a frequently 
occurring task in the inlegrational calculus. For that purpose 
one has great use of the following theorem due to Jordan (Cf. 
Appendix 2): 

Let CP (p) be a function integrable in the region Re [p) > a 
along any smooth curve, when | p | > c, and such that | CP (p) 
tends uniformly to zero when |p|-^cxD. Let C denote a half- 
circle with radius q and centre at (cc, 0) in the same region. 
Then 

Urn j (P (p)-c^^-dp — 0 (8) 

a 


for all negative values of t. 

Naturally a corresponding lemma subsists for positive values 
of t and an infinite half-circle convex to the left if (P (p) fulfils 
the same conditions in the region Re {p} < a. 

The equality (8) is easily seen to be vafid also in the case 
where the integration is extended over only a part of the half- 
circle C. 

Now, the general method of treating an integral of the type 

O’ + i CXD 

a — ioo • 


is as follows: 

According to a general proposition already mentioned for the 


class of functions / {t) in question, the function 


F{p) 

V 


is holo- 


morphical in the region Re {p) > a, and there it also fulfils the 
conditions for Jordan’s lemma. Another fundamental theorem, 
due to Cauchy, permits us to change the path L of integration 
into an infinite half-circle convex to the right. Then Jordan’s 
lemma states 1 to be zero for all negative values of t. Of course 
we could know this fact d priori^ but our reasoning has given us 
a sufficient condition for an integral I (t) to have this property. 

The calculation of I (t) when t>0 necessitates a more 
complete knowledge of the analytic properties of F (p) on the 
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left hand side of the line Re{p} — a. In most cases has 

P 

sucli properties that Jordan’s lemma is applicable also in this 
half-plane. Then we may complete the line of integration by 
an infinite half-circle convex to the left, thus getting a contour 
C, inside which the integrand function possesses certain sin- 
gularities as poles and branch-points; we assume the latter to 
be conveniently connected by cuts making the function uniform 
inside C. In accordance with Cauchy’s theorem the integration 
along L may be replaced by integrations around the poles (resi- 
due calculus) and along the cuts. In a following chapter this 
procedure shall be abundantly exemplified. 


2. The notion of integrator. Processes of differentiation 

and integration. 


In the integral theorem to each function / (t) another function 
F (p) of the comple^i: variable p is associated in a unique way. 
Obviously this function F (p) is characteristic of and can be 
used to represent / (<). Hence the unit-function H (t), e. g., is 
represented by F {p) — l. 

Now we are going to introduce a simplified notation for the 
integral expression (1), writing 


f(t)^F(p)-H(t) (9) 

Here the expression F (p) symbolizes a definite operation of 
integration that is to be performed on the function H (t). For 
that reason we may call F (p) an integrator. As / (t) is obtained 
by the application of F (p) to the simplest and most fundament- 
al function H {t), we also will call F (p) the basic integrator of 
/ (t). Of course, any function of class 1 possesses a uniquely 
defined basic integrator given by formula (2), but also other 
functions may be given by their basic integrators. 

More generally, we shall call a function F (p) of the complex 
variable p a basic integrator, when the corresponding integral 
(1) converges and remains stationary for all values of a suffi- 
ciently large. Hence, by definition, the path L of integration is 

F ip) 

prescribed to be situated on the right of all singularities of — — . 
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Later in this chapter we shall take under discussion the inter- 
pretation of an integrational expression in some cases where the 
integral does not exist. 

By the symbolic notation 

0{p)-f(i) 

which should be read: »integrator © (p) upon / {t)y>, we under- 
stand the function [<!> {p)-F {p)]-H (t), if it exists. 

From the definition of a basic integrator we immediately con- 
clude the validity of certain arithmetic rules. Hence we have e. g. 

O{p)-[k-f{t)\^k-0{p)-f(t), 

®(p) • [A it)+h (<)] = ® ip) - A it) + ^ ip) -hit), 

[0l{p) -f ®2(p)] ■f{t) = ®l(p) •/ (<)+ ®*(p) •/ it) 

[©1 (p) • ©a (p)] •/ (0 = ®i(p) • [®2(p) 7 it )] = ©2 (p) • [®i (p) 7 it)] 
[©i(p)-©2(p)-©2(p)]7(0= 

©l(p) • ©2(p) • [©s(p) 7 (01 = ®4P) • [®2(p) • ©»(p) 7 (<)] 


More precisely, we find the so-called distributive, commutative, 
and associative laws for processes of addition and multiplication 
to subsist in the calculus with basic integrators. 

When applying different integrators to a given function 

/ (t) one obtains an unlimited number of n^w functions, which 
generally belong to other classes than the original one. Obviously, 
when (D (p) is itself a basic integrator, the operation is always 
possible. When 0{p) is not a basic integrator (e. g. <2)(p)= a 
polynomial in p), (D{p)-F(p) may be a basic integrator or not, 
according to the case. 

The properties of the integrator function F{p) depend of course 
wholly upon the primary function f (t). To get some knowledge 
of that relationship we will perform integrations by parts in the 
definition formula (2). As a>a we find: 

“ |< lim \ i ^ 

IpI 


Urn 




=0 


and the first step of integration gives the result: 


F(p)=/(0)+ 

0 


( 14 ) 
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If there are exponentially limited derivatives of higher orders, 
the same procedure can be repeated one or more times. Thus, 
when j{t) is of class n at least, one finds the following expres- 
sion for F (p): 


F(p)=/(0)- 




p 






pit 


pn 


(15) 


0 


The integral appearing in the last term, which we may call 
the rest-term, is obviously uniformly limited in the domain 

Re {p) > ao > «. 


If / (t) belongs to the class of order infinity, n may be any 
positive integer, and the rest-term will approach zero when 
n->oo, if |p|> 1. Hence, for such a function the following 
expansion holds: 


F{p)==fiO)+fM+W+. 

^ r r 


. 4-^)4. . 


in inf. (16) 


Here the right hand side constitutes a power series in ~ that, 

as we know, converges for some finite value of p. Outside any 
circle, concentric ^vith and including the circle of convergence, 
the power series converges uniformly, and F {p) is there uni- 
formly limited. Thus we have the general theorem: 

The basic integrator of a, function / {t) belonging to the class 
of order infinity is regular at the point at infinity. 

The properties of F {p) just stated make it possible, for ^>0, 
to complete the path of integration L by an infinite half-circle, 
convex to the left, and then to perform term by term the com- 
plex integral: 


C+L 


no)+m+m+ 

p p^ 




^pt pi 

Now, the residue of the typical term /^^^(0) -“ is 
and thus we find ~ 


/(<)= 




( 17 ) 
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which is the Taylor-expaiisioii for / (<). Any function of the 
class of order infinity may be represented in this way by a series, 
uniformly convergent in any finite interval of positive values of t. 

As a special result we get from these considerations the im- 
portant practical integrator rule: 

(18) 

for the present proved only for positive integers n. 

In the preceding pages we have seen how a function / (<) in 
a convertibly unique way is represented by its basic integrator 
F{p), Now the question arises how the basic integrators of the 
derivatives f\t), if they exist, are related to F (p). The 

formula (14), or more generally (15), will give us some informa- 
tion in that respect. 

If the function f (t) is of class (n+1) at least, /^ (<) has a 
derivative, piecewise continuous and exponentially limited; hence 
the integral theorem is applicable to (t) and the basic integrator 
of this function will be 



From the formula (15) we obtain the expression desired for 
Fnip): 

Fn{p)=p^- 

or 

Fn {p)=p^-F {p) — [j (0)-p"+/'(0)-p"^^H (0)-p] , , (20a) 


F(p)-/(0)- 


m /"(O) 


/<n^i) (0)' 


• ( 20 ) 


The last constituent of the right hand side, which we may call 
the principal part of p^F{p), has a convergence-producing task; 
by subtracting it from p^-F{p) a convergent integral is obtained, 
which of course must be the case as this expression necessarily 
represents a basic integrator. 

From (20) results how a linear differential expression with 
constant coefficients can be algebraized by the passage to the 
corresponding integrators. 
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On account of (20) we can pronounce the following proposition: 

To apply to the function / (t) the derivative symbol means 

the same as multiplying its basic integrator F (p) by p, provided 
that the initial value / (0) is zero. To apply the integrative 

t 

symbol J" dt means, without restrictions, the same as multiplying 

0 

1 


the basic integrator by 


V 


The latter part of the proposition follows immediately from the 

f 

former one, if we note that cp{t)— f j (r) dr is a function with the 


properties: (p{0) = 0, 

Using our mode of writing, when /(0) = 0, we have the equal- 
ities 


I / > {pyH{t)=p-F{pyH {t)=-p-f it) 


lass of 


expressing a certain eciiiivalence between the symbol -ierivative 


symbol p. To gCt a deeper knowledge of that weined oper 
attack the problem of differentiation from another side 
from tlie integral 


a+i oo 




?)•/(<) for 
I is not a 
erge, we 
ntegral, 


a — ioo 

By formal differentiation under the sign of integratioiiommu- 

G+i oo always 


1 r 


nay be 


thus a factor p before F {p). But in ma ty cases the new inte- 
gral does not converge, and then the ^jrocess of differentiation 
cannot be performed in this way_ r'or this a uniform conver- 
gence of the new integral is req«*^ed. It is, however, possible to 
modify the original integral h. such a way that the procedure 
is correct, provided that F (p) satisfies certain conditions. Leav- 
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ing aside the question of the most general form of these conditions, 
we will make some assumptions relative to F (p) that are usually 
fulfilled in the practice. 

F(p) 

We assume to be holomorphical in the region to the 


right of the straight lines p=a+r-c^'/’ and p=a+r e“^'/’, where 


It > 7;> — , and such that 


|^(P) 

I V 


tends uniformly to zero when 


p->-oo in this region or in any region contained by it (Cf. fig. 2). 
Such a function or basic integrator is said, for brevity, to satisfy 
the condition A. 



Fig. 2. 


As before we co dude that the function f (t)=F {p)-H {t) is 
zero identically when ’ is negative. When t is positive we may, 
on account of the pro^ ’rties stated about F (p), complete the 
line L of integration by pa. ‘s of an infinite half-circle, convex • 
to the left, and with its centre at {a, 0), and then, without tra- 
versing any singularities, change the path of integration into two 
straight lines Li and Lt, as shown by the figure. The integral 
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along Li and Li we will call the modified integral. Now, this 
integral is easily proved to be uniformly convergent for all t^U 
(ro== an arbitrarily small positive number)*, and this is also the 
case with any integral of the type 

^.■f WFip)‘^-dp (« 

Li Ls 


where Q (p) denotes a function that in the region in question 
increases at most as a power of p when |p|-^oo. 

From this we conclude that, when an integrator F (p) satisfies 
the condition A, the corresponding primary function / {t) for 
positive* values of t has an unlimited number of derivatives, and 

that the symbol ^ before / {t) means the same as the factor p 

before F (p) in the modified integral and conversely. One should 
observe that the new functions p-F {p), p^‘F (p) etc. as a rule 
do not satisfy the condition A or have the character of basic 
integrators. 

Thus it proves possible to establish for the specified class of 
functions F (p) such an equivalence between the derivative 
d 

symbol ~ and the* integrator symbol p that a well defined oper- 
ational meaning can be given to the expression ®(p)*/(^) for 
positive values of U also in cases where ^(p)‘F {p) is not a 
basic integrator. When the' integral (1) does not converge, we 
will always interpret this expression by the modified integral, 
provided that F (p) and (Z>(p) fulfil the conditions stated. 

At this point, however, we must emphasize that the commu- 
tative law, which is valid for basic integrators, does not always 
subsist for this more general class of integrators. This may be 
exemplified in a simple case. The expression 

phH{t) 

means, in accordance with the definition, if <>0: 


ph-H(t)-- 


1 r , 

. / pi i 

mj ^ p 


27ti 

LiLi 


dp- 


d^ 

dt^ 27ti 


K fvAv- 

mj 


V 




LiLa 


* Cf. Appendix 4. 


B 
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Thus one has the right to make the factorial decomposition 
pi=pa.pj and then to replace p® hy the symbol^. Quite as 
well, however, one could try to proceed in another way, viz 

pl-H(t)=p^p^-ff(<)=pl-^-H(<)=0. 

Obviously this procedure is wrong as it violates the integral of 
definition. Hence the commutative law cannot be applied to the 
product (i). 

As a general rule all operations of differentiation ought to be 
performed last in the suit of operations occurring. Indeed, by 
differentiations certain constituents of the time function can dis- 
appear that would perhaps have some influence by another mode 
of proceeding. In any case the correctness of the procedure can 
be checked by an examination of the modified integral. 

Making use of our generalized class of integrators we will 
now evaluate (/) according to (20 a) 

/<n) (0- vV-l^(p)-[/(0)*p^+f(0)-p^-^+ • • • {0yp]}-H(t) (22) 

The modified integral* gives, when t>0, 

LiLn 

( 23 ) 

(0)-p"+/'(0)-p"'''-i (0)-p]~-dp. 

Now, for the path of integration in question the following 
general equality holds: 



LiLi 


and hence, if i > 0, 



LiLi 


* Fn (p) is supposed to satisfy the condition A, 
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The disappearance of the last part of the right hand side of 
(22) for positive values of t will of course also occur if we in- 
terpret the symbol p as “ • 

Instead of employing the modified integral to interpret an 
integrational expression F{p)-H{t) when F(p) is not a basic 
integrator, one can make use of the practical rule of simply 
suppressing the principal parts of the expressions, which are 
often easily isolated, and then retaining the original path of in- 
tegration L. In fact, that is what was done, when in (20a) the 
parenthesis [ ] was subtracted from p^-F{p). 

The procedure may be elucidated by a simple example. From 
the foregoing we have the rules 

.. (26) 

p -f~0) p ’■f-O) 


Suppose we have to interpret the expression: 


ap*+b-p^ 


(27) 


The integral of definition 

/ 


(J-fi CK> 

1 pip^+bp^ ^ 
27tiJ p 


dp 


(28) 


has no meaning in this form. By a simple process of division 
we obtain 


ap*+ bp^ 


ap^+bp—ao)^ 


p®+ 


.... (29) 


If here we suppress the principal part and introduce the 
remaining constituents in (28), we get two convergent integrals, 
which, in accordance with (26), can be at once identified with 

— a-w^-cos (^) and -~6w-sin cjMf (^). 

Hence the function in question is 

/(<) = ( — aw^-cos (ot — bco ' sin a)t)-H (^) 


(30) 
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A simpler proceeding is as follows: 




(0 ^ 


■H(t) 


b 

la*-T;;TCOS o)t-i ^];;;‘sin cot]'H (t) 

av (0 at ' 


(31) 


==(— a-f//^-cos (ot — 6*wsin (ot)‘H {t). 


At last we will determine the basic integrator of / (t). This 
can be done in two ways. On one hand we obtain the integrator 
wanted from (29) after suppressing the principal part. Hence 

ap^ + bp 
p^+co^ 




(32) 


On the other hand, according to (30), F (p) is the sum of the 
basic integrators of — aio‘‘- cos o)t-H {t) and — bco •sin cot- H (t). 
Thus, from (26) we get back the equality (32). 


3. A linear integral operator and its. properties of 
composition. 

Let / (t) be a function of class 1 at least. We start from the 
identity 

t 

0 

which may also be written 


/(<)=/(0)-H(t)+ 


-/no- 


H{t-c)-dr, 


(33) 


and obtain, taking t—t as a new variable of integration. 


/(t)=/(0)-H(t)+ 


J'nt-T)- 


H{i)dT (34) 


The last identity can be written in a symbolic manner 


f{t)=[f{0)+J'r{t-Tydv]XH{t) (35) 
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Accordingly, with any function / {t) we can associate an integral 
operator F: 

F^[f{0)+ J'f{t-T)-dr] (36) 

0 


and generally define a symbolic product : 


Fxcp{t)=f{0)-<p{t)+J'f{t-T)-(p{T)-dv (37) 


Remark. The term » operator » is here reserved for such symbols 
as indicate a process of direct differentiation or integration to be 
performed on a function. The denomination » integrator » is em- 
ployed for symbols denoting that a new function shall be derived 
from one already given by some well defined procedure of inte- 
gration of a more indirect nature. In many cases the same expression 
can simultaneously symbolize an operator and an integrator. 


The common arithni.etical laws are applicable to the functional 
operator F just introduced. From its linear nature we at once 
conclude the validity of the first three rules given for integrators 
on page 12. To ^get the rules for symbolic multiplication we 
must examine the composition properties of the operator. 

Let /i (t) and /2 (/) be two functions with the corresponding 
operators 


Fi=[/.(0)+ J'j\{t-x)di] and F,=[f^{0)+ 


J dc]. 


0 


0 


We shall have to compare the results in the two cases of 
composition 

FiXiFoXjit)) and F 2 X{FiXf{t)), 


where f{t) denotes any function. According to the definition we get 

t t 

?> (<)=[/> (0)+y/’i(<-r)-dr]x[/.(0)./(f)+y'A(i-r)-/(0-d'^] = 

0 0 

t t 

=/-(0)-A(0)-/(0+A(0)- y' Mt-ryf{T)-dv+U{0). J~ \'{t-T)-f{v)-dx 

0 0 
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t T 

+ {'^i)'dvi. 


The double integral can be transformed conformably to a method 
due to Lejeune-Dirichlet (Cf. fig. 3). 



t t 

= f{Tl)dTiJ'fi{t — T)-fi{T — ti)-dc 


Hence, with the notation 


(t-T)=y/'> {§-T)-d§, 


and using a symbolic mode of writing, we obtain 

i^iX(FaX /(<))== (40) 

t t 

[A(0) /2(0 )+/i( 0)- J'Mt-Tydx+UiQ)- J'i\ (<-r) dr+ 


t 



(/— r)*d'r]x/ (<). 


0 


0 
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It is easily seen that fn (Ut) is really a function of the 
argument {t—v); indeed, a simple change of variable of inte- 
gration proves that a function of the type 


V 

F {x, y)= j'u (^— ®)-d| 


X 


satisfies the equality jP(£C+a, ^+a)=jF (a?, «/) for any value of a. 
It is also simple to verify that such a function F y) is sym- 
metric in u and u, i. e. remains unaltered when u and v change 
places. 

Hence we have the identity 

t 


(41) 


t 


Now, if we treat the second composed operator F 2 X(.Fi x/(0) 
in a similar manner, we must, by symmetry, get the same result 
as in (40), if we let the indices 1 and 2 change places. But, 
on account of (41), (40) is invariant against such a procedure, 
and thus we may write: 

F, X (F^ XJ ( 0 ) ^F,X{FiXf ( 0 ) = (Fi XF2)Xf{t).... (42) 

Accordingly, two integral operators Fi and F 2 , which operate 
consecutively on a given function, their order being immaterial, 
may be replaced by a new^ operator, which can be regarded as 
obtained from the former by a kind of symbolic multiplication 
in agreement with the principle: 


FiXF.^F^XFi^ 
t t 


[/i (0)4- J"fi(t—TydT\ X [/a (0) 4 - J^f’i (i— ‘r)-dr]= . . 


( 43 ) 
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We have just found the commutative law to be valid for this 
product. The associative law of multiplication, expressed by the 
equalities 

Fi X F 2 X Fz=Fi X {F^ X JFs) =Fi X (Fs X Fi) = jPs X (F 2 X Fi ), . . (44) 

can easily be proved in a similar manner. In fact, repeating 
the procedure employed, we get: 

X F2 X jF's = 


t t 

[/i(0)+yA(<-r)-dT]x[/.{0)+y/’.(<-T).dr]X[/s(0) + 


+ / A {t-r).dr\={U (0) /. (0)-/3 (0)+ 


(45) 


t t 

+A (0)-A (O)-y' '3 (t-r)-dT+U (0)-/i(0)- Z' f',(<-r).dr+ 

0 0 
t t 

+A{0)-/a{0)- y^A(<-T)-dr+/.(0)- J'rn{t-T)dr+ 

0 0 
t t t 

+A(0)- J'fn(t-T)dr+f,{0)- J'j\ 2 {t-T)dT+ J'i\.^{t-r)-drl 


The meaning of /Vs (<— t) (r, 5==1, 2, 3, r ^ 5 ) follows from (41). 
The function / 12 Z (^— ^) can be defined e. g. by the integral 

t 


/W(<-A = f 


(46) 


Using the method of Lejeune-Dirichlet one can at once verify 
that 

/'l23 = / 02)3 = /V23) = / 2(13)= etc (47) 


Thus the indices may be arbitrarily permuted and associated. 

It is often more covenient to write composed integrators under 
the original form (36). Then, instead of (43) we obtain 



25 


where 


t 

i^iXF.=[A(0)-M0)+ j'K^{t-T) Ax\ (48) 


(®)=/i ( 0)- a {x) +U (0 )-/’i (a;)+/'i. (a;)=F2i(a:) (49) 

We call Kn{t—‘T) the integral kernel of the operator F\XF<i- 
To illustrate the general rule of forming the integral kernel we 
represent the expression of the product of four factors: 


Fxxf.xFsxf.=[A(o)-mo)-/3(o)-mo)+ Jk 1234 {t — r)dr] (50) 

0 

Ki2U (1^) = 

+W74(0)-/i(0).M0)+A(a^)./.(0).M0)^ 

+ /-(^)-/3(0)-/4(0) + A3(^)-/2(0)./4(0)+^^ 
+/'2a(a^)-/i(0)-/4(0)+/24(a^)-/i(0)-/3(0)^ 

+ A23(a?)-/4(0) + /'234(a:)7l(0) + f34l(aj)-/o(0) + /'412(a:)-/3(0) + 

H~/ 1234 (x). 


As wee see thi:; expression will be considerably simplified if 
some of the initial values /,(0) are zero. When they all vanish, 
we simply get: 

t 


Fi X F'i X Fs X Fa — [ / / 1234 {t — T ) • dr] 


(50 a) 


0 


Of special importance is the case where all the operational 
factors are equal, that is when the product is a power of the 
symbol F. Then we conveniently use simple indices for the 
functions f and K. The original function being / (t), and 
accordingly 

t 

F=[/(0)+y/'(i-r)-dr], 

0 

we obtain the following general expression for an integer power 
of F: 
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where 


K 


ft 

F^=\i (0)»+ J' sr„ (<- t)dTi 

0 
n 

ysst 1 

t 

i\{t-T)= J' ^'n-.{t-§)-f{§-r)-dS-, i\{x) = i\x). 


The proof is readily carried out from n to n+\ in accordance 
with the mode of proceeding just employed in the general case. 


4*. The equivalence between integrators and integral 
operators. The inverse operations. 

According to the foregoing we can associate with any function 
/ {t) two operational symbols, the integrator Fijp) and the integral 
operator F, which, if applied to other functions, give rise to 
new ones. The definitions were: 

oo t 

F {p)=f (0)+ y' {r) er-P^dT; F=[f (0)+ j''r {t-r)-dA. • (52) 
0 0 

Operating upon the unit-function, H {t\ both give back the 
original function / (t). Operating on an arbitrary function (p {t), 
they will both give the same new function, that is: 

F(pyq){t)=FX(p(t) (53) 

To prove this identity we will make use of an auxiliary 
theorem, the correctness of which immediately follows from the 
definition of an integrator, viz 

(T+too 

H ® (p) “ • dp=[0 {p) H («)] (54) 

a — ioo 

At first we transform the right hand side of (53), utilizing 
<p (t)'s prop erty of vanishing for negative values of t 

* Cf. Appendix 7. 
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F X9)(0=/(0)-95(«)+ J'f{r)-<p{t—TydT 
0 

oo 

=f{0)-(p{t)+ J'f(t)(p{t—r)-dT 


=/ (0)-9’ (<) +y /' (p)H(<)]<-r-dT. 

0 

Hence, on account of (54) and the integral theorem (Cf. Appendix 1 ), 

oo 

FX(p(t)=rf {O)-0{p)-H{t)+ f (p).H {fy dx 


OO 

=[/(o)+y/' 


{t)- • dry 0{p) H {<), 


and at last, according to (52): 

FX(p{t)=F ip)- 0(p)-H {t)=F {p)-<p {t). 


We also know the symbols F (p) and F to obey the common 
arithmetic laws of addition and multiplication, and therefore we 
can regard them as equivalent when simultaneously appearing 
in expressions obtained by these operations of calculus. This 
fact is of great practical importance for the following reason: 

The employment of integrators for solving differential equations 
will present the unknown function under the form 

/ (t, .)=F {p; x,y,z.. .)-(p {t) (55) 

To find the explicit function of time one can of course always 
fall back upon the definition-formula (1) and try to evaluate 
the complex integral. In many cases, however, the function 
F ip; 0 ?, y, 2 . . .) can be broken up into terms and factors in 
different ways, and these simpler constituents are often of similar 
types from one case to another. Thus, having once derived the 
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integral operators, corresponding to a number of such typical 
integrators, one can, on account of the equivalence just mentioned, 
make use of these results and in this way often simplify the 


interpretation of composed integrational expressions. 

Assume e. g. we want to find x (<) when 

X (t)=Fi {p)-Fi {pyFt (p) F4 {p)-<p it) (56) 

and that we have already obtained the equalities 

Fi (p)-H {t)=F^xH(tl F, [pYH (t)=FsXH{t). 

Then we rewrite x (t) in the form 

x(t)=F2{pyF,{pyFiXF,X(p(t), (57) 


and have only to evaluate the complex integral 

a+i<x> 

a — ioo 

This brings us the operator 

U=[u (0 )+ i (t-r) drl 
0 

Hence, the final computation of x (t), ^ 

x{t) = UxFiXF3X(p(t), (58) 

according to the rules of composition of integral operators, is 
lead back to pure quadratures. 

The method of calculus described above is not always the 
most convenient from a practical point of view. It is often 
better to retain the form (55) of the integrational expression and 
sometimes even to introduce 0(p)-H{t) instead of (p{t). This is 
especially the case, when F (p)'0 (p) is of such a form (e. g. 
is a rational function) that one of the so-called expansion theorems 
can be used. (See later in Chapter 7.) 

The relationship between a basic integrator F{p) and the 
corresponding operator F can be elucidated from another side 
by a study of the inverse operations. 

Proposed that (p (t) is any given function with the basic integrator 
0 (p), on applying F (p) to tp (t) we obtain a new function, say yj (t) 

xp(t)=F{p)-(p(t), 


(59) 
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the certainly existing basic integrator of which is 

>F(p)=F(p) ® (p) (60) 

The function y> {t) can also be defined by means of the integral 
operator F: 

ip {t)=FXqjit)=f {()) (p{t)+ f'{t—t)-q){T)-dt (61) 

0 

Our examination now refers to the opposite or inverse operations, 
which bring us back from ip{t) to (p {t). From (60) we get 

1 




(62) 


Therefore 


= U {p)-ip(t) (63) 

Provided that U (p) is itself a basic integrator, employing its 
corresponding operator U, we may write (p (t) under the form 

• t 

(f {t)==UXy){t)=ru{0)-y){t)-{- j* (64) 

0 

Now, the equations (64) and (61) being reversed to one another, 
we ought to establish a functional relationship between fit—r) 
and u{t--T) by inserting the expression (61) of y)(t) into (64). 
By this means we get the \dentity: 

t 

(p (<)=« (0) •/ ( 0 )- 7 i (<)+« {^)- J'f{t—r) (p (t) dr+ 

0 

t T 

+ J' dT-[j {0) (f{T)+ J" '’{r—T,)-q{Ti)-dTi]. .. (65) 
0 0 

The double integral may be transformed in a well-known way, 
and thus we obtain: 

[l-/(0) tt(0)] (p(<)= (66) 

t t 

J p (T) dT-[u (oyf(t-T) +/ (0)-«'(<-r) + y' ■f(i-r)-dil 
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Now, if this homogeneous integral equation shall be satisfied 
for arbitrary functions (p (i), the following two conditions must 
be fulfilled: 


/(0)-«(0) = l; 


u 


I 


..(67) 


These can be combined to one, if we introduce the notations 


( 68 ) 

Then, instead of (67) we get the only condition 


/*(^— r)+^*(f 




'‘*'(^— ^) * r) • = 0. 


(69) 


We will call U the inverse integral operator to F. As for any 
function (p{t) we have q) {t)=:UXFX(p {t), the relation I7xJP=l 
subsists, and instead of U we may conveniently use the notation 

(70) 

According to (67) a necessary condition for the operator F to 
possess an inverse operator is that / (0) =i= 0. Then also u (0) =f= 0, 
and. the functional relation (69) subsists between the functions 
u* (t^r) and /*(<— t). Two functions, linked together by this 
relationship are said to be reciprocal to each other. The equation 
(69) is another mode of writing the relation F(p)*[7(p) = l 
between the basic integrators or the relation F X [7 = 1 between 
the integral operators. It constitutes an integral equation of 
VolterroLS type. If we regard /*(t/) as the given function and 
{y) as the function to be determined, we can e. g. rewrite the 
equation under the standard form 

y 

/7*(2/-~iJ?)-'W*(a?) da?—— /*(i/) (69 a) 


and then obtain the unique solution by forming the resolvent 
according to a wellknown general method. 

More interest is offered by the expression for the inverse 
operator F^S which is easily derived in a similar way. One 
will find, by solving the Volterra equation (61), that 
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F 


t 

-‘=[/(0)->+yF 


(t-T)dT\ (71) 


where 




f.{x) 

/(O) 


y-rl 


.... (71a) 




The functions fy{x) are obtained from (51). 


As a simple example of the inverse operations we take the 
integrational equation: 

(“) 


Here the operator is 


+ co^ 


-F=[l— / sin w (<— r)*dT]. 


Regarded as the unknown function, {t) is determined by the 
integral equation 


(p{t)—oj- j sin CO {t—T)*(f {r)’dT= yj (t) 




On the other side 9 ; (t) can also be obtained from the equation 
(a). Hence 


CO' 


( 0 - 

Here the right hand side can at once be interpreted, and we get 


t T 

(p{t)-=xp{t)+CipJ' /'(ri)-dri 


or 


— T)-y>{T)-dr. 


iy) 
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Therefore, the inverse integral operator of F will be 

t 


F —^ = [ 1 + — -r) • dr]. 


It is easy to check that FXF~^=1, e. g. by inserting the 
expression for y}(t) according to (/?) into (r)] a pure identity is 
then obtained. 


Our considerations furnish us with a general method of solving 
a Volterra equation of the special type: 

t 

q){t)+ f'K{t—T)q){t)-dT = ip{t) (6) 


This becomes identical with (61) if we put /(0) = 1 and 
K{x)=f{x). On account of the first of the formulas (52) the 
corresponding basic integrator is 


ii^(p)=l+ K{Tye-^HT. 


(^) 


If this integral is comparatively simply evaluated*, it is in many 
cases practical to present the solution of the integral equation (()) 
in the form of the integrational expression 

{»?) 


In our example we had K {t) ——(j) siiiwr; hence 

oo 

F{p) = l — w / iin o)T-er-P^dT—-~, — r, 
J p^+O)^ 

0 

and the equation {rj) gets a very simple form. 


* It may happen that K {^) can be approximated by a suitable aggregate of 
polynomials and exponentials. 
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5. A table of some important integrational expressions and 
rules of calculus. Applications of the general methods of 
the foregoing chapters. 


In the table below we present a collection of some simple 
integrational formulas that are of frequent use in the practice. 
This collection will be subsequently completed in the following. 


1 . 


2 . 




3. 


4. 

5. 


_? --Hit ) — 


(p-«y 


1 


H{t) 


p—a 


6 . 


•p^+ (id ' ' ' ' 


n (p+/?)-sin 9)4 -w • cosf/;' ^ rr/.x 

7. >-H{t)=e 


CP {p+a)-j {t) 


8 a. (P (p) /(<)=c"^-(P(p+a)-c”®^*/(0 

9. c«^-/(0=/{^+a) 

10. If F{p)-H{t)=f{t), then F(Ap)-H X>0 


The formulas nrs 8 and 9 are sometimes called the ^^shifting'"- 
rules, nr 10 the **uniformiUj'-rule. 

All these formulas are easily proved by means of the integral 
of definition (1). Some of them, viz nrs 1, 2, 5 and 6 have 
already been deduced. As instructive exercises upon the general 
methods of the preceding chapters we are now going to verify 
and discuss several of them in different ways. 


c 
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At first we compute the integral of formula 7 by residue 
calculus, after completing, for i>0, the line of integration L by 
an infinite half-circle convex to the left. The poles are p=—fi±i<o, 
and integrating around these we obtain 



iwsin cos (p 


, 0 — 


i 


= — •l?^{(fosin 9 o — i(o cos(/;) -)-(/;). 


This rule can also be deduced from the formulas 5, 6, 8 and 
2. First we get: 

sin {t) = cos 99-sin dd-H (i^) + sin (/>*cos (ot-H (t) 

_ p-sinf/)+fr)-cos9) 

=p — H{t) 

Then, according to 8, we can make the transformation 
sin {o)t+(p)^H{t)^-0(p)-H{t) = a){pyef^ie-^^>H(^^ 

Hence 

sin H {t)^(D p-\- (i)' 

Inserting here the expression for d) (p+p) we get formula 7. 


The rule nr 8 can be proved e. g. in the following way: 

(7 + 1 oo 

®(P+«)/(0 = ®(P+«)-^(p)ff y' n^+ayF{A)-^^dA. 

a—i oo 

By the change of variable X+a—p we get 

( 7+100 

^(p+a)-f{t)=er^-~^J 0{p) F{p—a) -^-dp 

a — ioo 

keeping the limits of integration unaltered. This is permitted as 
the value of the integral is invariant against a displacement of 
L to the right or to the left, provided that L is always situated 



sufficiently far to the right. According to the rules for integra- 
tors, ®(p) can he put outside the sign of integration, and thus 
we obtain 

<r+ioo 

ff—ioo 

again denoting by A the variable of integration. At last we 
introduce p ==/-—«, getting 

a+ioo 

1 r 

0{p+a)i F{p)—dp, 

a—ioa 

that is 

O (p + a) • / (<) = • (D (p) • • / (f) , 

which is the result desired. 

The function j{t) being arbitrarj% we only need to replace it 
by in order to pass from the formula 8 to the for- 

mula 8 a. 


It may be of some interest to show how the rules nrs 1 and 
8 can be deduced from the integral operators of the original 
basic integrators. • 

In the former F{p) = —, and the corresponding integral opera- 

V 

lor, which we may denote hy /, accordingly will be: 


t 



0 


Therefore 


P 





T2 Ti 

dri /^(r)-dT. 



We can transform this repeated integral of order n into a 
single quadrature by making use of the integral kernel corre- 
sponding to Obviously, this kernel may be at once obtained 
from the first rule, but here we prefer to deduce it from the 
general formula (51). Putting f(x) = l we get by this 
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Mt—r) 


t 


fn (t—r)-- 




\n — 1 


Hence, as /(0)=0, the integral kernel becomes 


and accordingly 


Kn{^x)^fn{x)=^-, 




Thus we have got the formula 

t 

t—rT-^-q'{T)dT. 

0 


(73) 


If here (p{i) is replaced by H{i), the first, rule of the table is 
obtained. 


In order to derive the third formula in a similar manner we 
will start from the integrational expression 

(p—aY 

From rule nr 4 follows that 

and consequently, as /(0)“0, 

t t 

F»= Kn (t-T) = I y' 


Further is 


fi (t—r) fn (t — t) = I •(<— 
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Therefore 

t 

(p ( 0 90 {i) = (tydt. 

0 


If here we put fp{t)=H [t) and take u=t — r as variable of in- 
tegration, we get for i > 0 


(p — Of)” 


H{t)- 


\n--l 




Now we only have to add the factor p=-v: to both the mem- 

dt 

hers to recover the rule in question. 

Of course the employment of residue calculus will bring us 
the result in the simplest way. Indeed, from the expansion 


^pt _ ^at , ^iip — a) _ ^ut 


1 + |Y-(P a)+Y2^p—aY+- 


] 


— 1 

we find at once p 0 “^ to be the residue of the pole a, 
I n — 1 ^ 


At last we will give some examples of different ways of treat- 
ing a simply composed integrational expression in order to elu- 
cidate further the relationship between integrators and integral 
operators. 

If a sinusoidal electromotive force is suddenly applied to an 
inductance coil, a current x {t) will arise according to the 
equation 




0)-p^ 


(P + P) 


H(t) 


(74) 


To obtain the explicit time function one can e. g. proceed as 
follows : 
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A. Complex integration. 

For positive values of t we have 

0 

The poles are p = — /? and p — ±ia). Thus, by residue calculus, 

® W="| . -e^A 1 

\{^+ico)-2ico jj 

and accordingly 

x{t)=:^--^^^-{--p^e-^^+^-coscot+u)%\ncot)^II{t) (75) 

B. Partial fraction decomposition. 

By well-known methods one finds 

^ ^ I R V 

(p+)3) (p^+w^j ^ V~\’P ^ p^+w^ p^+w^y 

Then x (t) is readily obtained from the rules 2, 5 and 6 in the table. 

This manner of proceeding illustrates the use of the so-called 
expansion theorem, which, if applicable, will generally offer the 
simplest way of interpreting an integrational expression. The 
expansion theorem shall be completely treated in a following chapter. 

No^e. As a variant of this method can be mentioned tlye expansion of the inte- 

1 

grator into a power-series on One gets the Taylor expansion of the time 
function by applying rule nr 1 to each term of this series. 


We write 


C. Factorial decomposition. 


* W=^ (I) 

According to rules 2 and 6 we have 

(p)-fr(f)=e-/^^-JT {t), F^ipyH (<) = sin cot H (t), 
and consequently the corresponding operators are 

t t 

F\=[l — e(^^dr\ • * 

0 0 

Hence, on account of (48) and (49) 

® (t)=F,XF,XH (<)=[/i (0) /2 (0) + j' Kyit-x) dx]xH{t), 


(76) 

(77) 
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and the integral kernel is 


ifia(ir)==~-^^2-(r/j2-cos (78) 

Therefore we obtain the following expression for x{tY 

t 

J' [f/j2-cos wt^-“Oj^*sin 
0 

which is easily identified with the former (75). 


D. Mixed method. 


Employing the same factorial decomposition as before, we can 
alternatively present x (t) under the forms 

x(t)=FiXF 2 (p)-H(t) and x (t)=F 2 XFi (p)^H (t). 
Conformably to (76) and (77) these symbolic equations have the 
meaning: 


X (^)=(sin r/>r-dr)*JT {t) 


(79) 


0 


t 

X wT-dr)-// (<). 


(80) 


Of course, also these two expressions are but different ways of 
writing (75). * 


6. Application of integrators and operators to the solution 
of differential equations. 

As already emphazised in Chapter 1, the practical applica- 
bility of the notion of integrator to the solution of differential 
equations requires these to be linear and to have constant 
coefficients with respect to the variable of time. In fact, in this 
case the introduction of integrators makes possible the more or 
less simplifying procedure of algebraizing. 

One can regard the integrational calculus as a generalization 
of a well-known method using complex vectors to find in a 
simple and lucid way those solutions of differential equations of 



the kind mentioned which correspond to stationary sinusoidal 
disturbances. In the alternating-current theory this method of 
calculus is usually called the ''symbolic method"'; there it has got 
an immense practical importance, making it possible to apply 
formally the notions and laws of the direct-current theory to 
steady-state sinusoidal processes. By the integrational calculus 
this will be possible also to transient phenomena. 

In the symbolic method every lime function is put under the 
form where denotes the representative vector of 

the magnitude in question. The exponential factor can always 

be removed from the equations, the symbol ^ being replaced by 

ico and the symbol ^ dt by 1/ioj. (The initial conditions are 

left out of consideration). In this way one gets to solve an equa- 
tion or a system of equations no longer containing the variable 
of time. At last one finds the explicit form of the quantity desired, 
on multiplying the solution of this system by and taking e. g. 
the imaginary part of the complex number thus obtained. 

In the integrational calculus the simple sinusoidal function 
has its counterpart in a class of functions which are charac- 
terized, as we know, by the properties of being all furnished 
with the factor H (t), expressing their identicgil disappearing for 
negative values of t, and further of possessing a basic integrator 
defining them uniquely. This latter symbol F (p) can be regarded 
as corresponding to the representative vector F in the symbolic 
method, just as the unit- function H\t) can be said to correspond 
to the exponential factor From a methodical point of view 
we can say that the above mentioned procedure of obtaining the 
explicit time function in the symbolic method is quite analogous 
to the forming of the product F{p)-H{t) and the evaluating of 
the integral of definition (1) in the integrational method. 

We recapitulate the three main properties upon which the 
applicability of integrators in solving diCFerential and integro- 
differential equations is based, viz 

1. The law of uniqueness. 

Two functions with identical basic integrators are themselves 
identical and conversely. 
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2. The law of differentiation. 

The basic integrator of the nth derivative /<">(<) of a function 
f{t) with the basic integrator F(p) is: 

p"F(p)-[/(0)-p»+/'(0)-p«-*+ • • • • +/<’‘-*>(0)-p]. 


3. The law of integration. 

The basic integrator of the function 

0 

being the basic integrator of f(t). 


f{r) dr is ^ Fip), F{p) 


As a first example of the practical usefulness of integrational 
methods we choose the common ecjuation of oscillatory systems 
with one degree of freedom, 

(81) 

and assume the initial conditions 

(llT 

x{+0)=u and (+0 )==t;. 


Writing x{t)—X{p)'H{t) we obtain, on account of the law of 
differentiation, 


2^-{pX-pu)H(t), 

^^,=(p-X-p-u—pv) H (t) 


(82) 


Now, by introducing j {t)—F {;p)'JI {t) and employing the law 
of uniqueness, we get the equality 

(P“+2/3’p+o/“)-.Y {p) = F (;p)+{p-+2pp)u+pv. 


Therefore the required solution of (81) will be 
x{t)^Xi{t)+X2.{t\ 

where 


Xi{t)- 


iL{p)_ 

’p''*+2//p-t-w- 


H(<) = 


1 

p“+2/jp+(ir 


■m 


X.,{t) 


Mp°+(2/?M+r )p 

p-^+2/lp+or 


.... (83) 
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If F{'p) is previously known and of a simple constitution, as 
e. g. in the case: 

/(^)=^o-sin F{p)=eo--^—* 


the first of the two forms of is the most convenient for the 
following treatment. 

We may easily interpret the integrational expressions obtained 
and write down the explicit time functions if at first we de* 
termine the auxiliary function 




_P 

'p^+2(ip+( 




P __ 


■H(t) 


.... (84) 


( 01 ^ = 10 “ 


According to rule nr 7 in the table 

1 


a{t)= •c~/^^-sin o)it-H{t). 

(Oi 


(85) 


Introducing for a moinenl 

t 


.. ( 86 ) 


we can write the corresponding integral operator 

t t 

B = j^ft(0)+ yi»'(i-r) dr| = wi(i-r) drj. (87) 

0 0 

Hence we get 

t 

<2?i — r)*/(r)*dr 

(OlJ 


or 


t 

0 


( 88 ) 
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Further, as for />0 the factor p means the same as the sym- 
bol x^{t) can be obtained according to the equation 

a ?2 {< ) = — + 2)3 • j • — * • sin (OiU 


from which 


072 {t) = 


V + fjU . 




sin (Oit+u co^ (oit (<).... (89) 


111 studying transient phenomena in electrical circuits built 
up by a finite number of lumped resistances, capacities, and in- 
ductances one will arrive at systems of integro-differential equa- 
tions of the type 

n 



(u = l,2---n) 

Indeed, the network can always be divided into a number of 
meshes, coupled to one another by mutual conductivity, capacity 
and inductivity, a;id to each circuit Kirchhoff’s law may be 
applied. In (90) denotes the current and the electromotive 
force in the M mesh. Lar and are coefficients of resis- 

tance, inductance and capacitance respectively. We will assume 
the system not to be in equilibrium but without initial charges 
on the capacities at the moment when the forces attack it.* 
Introducing the integrators of the currents and forces, one gets 
the following linear algebraical system for the unknown basic 
integrators ; 

n n 

yZ,„,-It,{p)=^E,,(p)+S^p-L,„ri,.(0), (91) 

)•=■ 1 )--= I 

(,u = l, 2- • n) 

where 

Ztit>=p-Luy+Ru>’'i-—-yT~ (92) 

' ‘ P'^ uy 

* The effect of initial charges must be the same as that of constant electro- 
motive forces. 



44 


Denoting by D(p) the determinant |Z^y| and by Dfip{p) the 
minor, corresponding to Zfjy, we can write 

2.(p)=iv(p)+n(p). 1 


,14 = 1 


.... ( 93 ) 


n n 

Jy" ( ^^."*'** (P)- 

k=l 


(*/=!, 2--* n) 

Thus the integrators appearing are all rational functions of p 
if, as usual, the Eu(p) are rational, and obviously this must 
always be the case in problems leading to systems of ordinary 
linear differential equations with constant coefficients. As we 
shall devote a following chapter to the treatment of this prac- 
tically most important class of integrators, we may leave aside 
for the present a nearer analysis of the results gained. 


In applications to partial differential equations one meets 
integrational expressions of various kinds and, as a rule, of a 
more complicated nature. In order to illustrate how a pro- 
blem involving initial as well as boundary conditions can be 
treated by integrational and operational methods we are now 
going to make some researches on the physically and technically 
very important wave equation. In this chapter, however, we 
will only discuss the one-dimensional case, but we will start 
from the most general form of this equation, represented by the 
universal equation of telegraphy. 

Thus we take under consideration a homogeneous transmission 
line, the primary constants of which may be : r = the resistance, 
a = the leakage, I = the inductance and c == the capacitance, all 
per unit length The distance between the branches may be d. 
We assume the line to be exposed to an external magnetic field 
with the component B {x, t) of induction normal to the plane of the 
line and to an external electric field with the component E (a?, t) 
in this plane, normal to the conductors. 
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V (x,t) may be the voltage between the branches at a point 
on the distance x from one end at the time t, and i («r, t) the 
corresponding current. The total length of the line may be s. 
Then we have the following system of differential equations: 


flv j di , j OB 

— htt- 

dx dt dt 


di 

dx 


= a'V+c 


dv 

dt~ 


-cd 


dt 


(94) 


According to the fundamental rules we get from these equations 
for the basic integrators V = V{x; p), /=/ {x; p), and B=B{x; p), 
E=E {x; p) a system of ordinary differential equations, viz 


—-J^=r I+l p-[I—i (x, 0)]+d p \B—B {x, 0)] 
—~—af'+c-p[r—v{x,0)]—cdp[E—E{x,0)] 


.. (95) 


Note, To simplify our mode of writing we occasionally let capital letters 
mean basic integrators and not, as usual, integral operators. 


We make use of the abbrevialions 


O{x;p)-=dpB{x;p)—p[li{x,0)+dYi{x,0)l 1 

W(x; p)— '-d c p E {x; p)—c p[v{x, 0)—d ¥^(x, 0)]) 

thus getting (95) under the more lucid form 


df' 

dx 


=='{r+pl)I+(D, 


-|=(a+pc).F+.p. 

Hence, eliminating one of the unknown (|uantities 
d^r 1 

where 


6^={r+pl) (a+pc) 

—Fi {x; p)=(r+pZ)- ’P(a;; p)—0'x{x; p),| 
— Ei {x; p) — (a+pc) • (D{x; p)— ^P'x {x; p). j 


(97) 


(98) 

(99) 

( 100 ) 



46 


Thus we have reduced the original problem to that of solving 
an equation of the form: 

^-d=(p).r=-F(a;;p) (101) 

under some given boundary conditions. If Z\p) and Z'\p) denote 
the impedances of certain energy-sources at the beginning and 
the end of the line with the electromotive forces d (t)=E' {p)-H (t) 
and e\t)=E’'{p)-H {t) respectively, we gel the boundary conditions 

F(0; p)=E\p)-Z\p)-I (0.- p), 1 

V(s; p)=-£"(p)+Z''(p)./ {s; p),t 

assuming, as usual in the practice, the initial values of the 
currents and charges of the end systems to be zero. 

Now, it is not difficult to prove that it is always possible to 
form the solution of the general transmission-line problem by 
superposition of the solutions of two more special problems, viz 

A. Homogeneous equations of the type (101) with jP' (a?; p) = 0 
and inhomogeneous boundary conditions: 

where: 

E'*(p)=E(p)+^^®(0;p) 

B. Inhomogeneous equations of the same type and homo- 
geneous boundary conditions like those above with 

E ^—E :i: = 0. 

The solution of the problem first mentioned is easily effected 
since the two constants of integration can always be determined 
in such a way as to satisfy (103). One finds 

V(x;p) = F\x;p)+r\x;p), 


(103) 


(103 a) 


( 104 ) 
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where 


shS{s—x)-\r 2 ^(s—x) 

r¥: P) = , ^z'z"\— r ,^+z“T7 ' ®’* 


(104 a) 


1 {x, p) . 2'Z” 


£a 


, Z'+Z” . v 

1 + 0 5 +^ ‘Ch 0$ 


E".(p) 


y__i/r + pl 

" V a+»c' 


In order to gel a general solution of the latter problem we 
conveniently introduce a so-called Greens function K {x, p), 

according to the definition: 

1 . K {x, p) is continuous in the interval (0, s), where, only 
the point x==^$ excepted, it satislies the homogeneous equation 

2. K (x, V?; p) satislies the homogeneous boundary conditions 

(103) with -0. . 

3. The first derivative with respect to x, K'x, is continuous 
over the whole interval, the point x—§ excepted, where we shall 
have: 


(IK{x,^:p) j 

()X Jl— 0 


Now, the Green function being found in one way or another, 
we obtain the solution of the inhomogeneous equation (101) b\" 
a simple quadrature according to the formula 


V(x;p)= / K{x,S;pyF{^;pydi, 
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The general expression of the function K {x, §; p) proves to be: 

vj ■ d{s—i)+Y-ch d (s—i) J 


KixJ;p) = 


Z' 

shdx-\-^-chdx] 
Z’Z' 

1 + 


Z'+Z" 


h d 2 
when X < ^ 


ch b s 


shd{s—x)+^„chd(s—x) \-\shd^+^-chd^ 

^ _ 


1 + 


7 ^ 

7 ? 


sh b s- 


Z+Z" 


chib s 


(107) 


when 


Finally, to obtain the explicit time functions we have to interpret 
the integrational expressions which we obtain when applying 
V{x;p), given by (104) and (106), to H (t). Taking under consi- 
deration only the most complicated case (106), we are thus faced 
with an expression of the form (Cf. Appendix 8): 


v{x,t} = V(x,p)-II {t)=^ / K{x,i;p)-F{s;p)-H (108) 


Here, of course, the expression under the sign of integration 
can be treated in the bulk according to the general integrational 
rules, and sometimes this is also the simplest way to arrive at 
V (x, t). However, as F {§; p) corresponds to a disturbance of 
external nature, while K{x,§;p) is an integrator that onlj" depends 
upon and to a certain degree characterizes the system in ([uestion, 
it is, principally at least, more adequate to treat each factor 
separately, i. e. to make use of the associated integral operators. 
Thus we may first determine, e. g. directl)^ by the integral of 
definition (1), the function 

k(x,^;t)=K(x,§;pyH(t) (109) 

which brings us the operator 

K{x,§)=[k{x,i;0)+ ^k't{x, t—T)-dr] (110) 

0 

Then, if F (^; p) is a basic integrator, we can derive its primary 
function 


f(i,t)==F{§;pyH{t). 
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This being done, we obtain the solution required under the form 


v{x,t) = 




( 111 ) 


or, more completely written 


S 8 1 


In the case where F{^; p) is not a basic integrator, the method 
described can always be applied to those terms or groups of terms 
in jP (if; p) which have this property (Cf. formulas (100) and (96)); 
as to the other terms one must, as we know, replace the factor 

p by “ and postpone to the last the effectuating of occurring 

CLl 

differentiations. 

Accordingly, the main problem will be the determination of 
the primary function k{x,§;t) of K{x,§;p), The universal ex- 
pression (107) tells, us how complicated this task can be. (In the 
practice the impedances Z* {p) and Z"(p) are usually rational 
functions of p). Frequently, however, one will have to deal with 
more special cases, which are more easily accessible to the ana- 
lysis. Thus, in a great number of problems the expression for 
6 (p) is reduced to a simple form, which will cause considerable 
facilitation. In this respect we may distinguish the following 
typical cases: 

1. The universal line. 

1 

Z=i2-(p+Mi)2-(p+^M " 

Here the notations mean: 

w = = R=\/^ (^^ 3 ) 

V/Zc ^ c \ c 

D 
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2. Line without energy losses. 

Here we have r=a=0, hence ^<i=t/2=0 and 

Z=R (114) 

3. Line without distortion. 

Such a transmission-line is characterized by the property 

r a 
L c 

Accordingly, 

d=--(p+M). Z=R (115) 

u 


4. Cable-type line. 

So we call a line the inductance I and leakage a of which 
can be neglected. In this case we obtain 

d=}/r'c-}/p,Z=\/--^ ( 116 ) 

▼ c Vp 

In the practice, of course, these last cases will only be ap- 
proximatively realized. 


Also the Green function can be considerably simplified by 
special boundary conditions. We will present some practically 
important examples of this. 


A. Line balanced at both ends. 

This means that Z'—Z”=Z. Accordingly, by (107): 

x<§ 

K{x,§; x>§ 


( 117 ) 
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B. Line short-circuited at both ends. 

Then Z'=Z" = 0 and we get from (107) 

.. , shdx- shd(s—^) 

K P) = 


K(x,§;p) = 


6 sh ds 
shd{s—x) sh 
d-shds 




(118) 


C. Line isolated at both ends. 

Then Z'=Z"=oo, and we get the same expression for K{x,i;p) 
as in B, except that in the numerator the sign sh should be re- 
placed by ch. 


At last we consider the special case of formula (104) where 
JS";}: = 0 and the line is balanced at the end or (which will prove 
to be principally the same thing) when it is of an infinite length. 
Then the integrator of the voltage function will be of the form : 


Hence : 


V (x; p) = 


V {x,t) — 


p—Sx 

- yrE{p) 

(119) 

p — ^x 

yre{t) 

>+l 

(119a) 


Many of the above mentioned special cases are mathematically 
equivalent to some classical problems in other domains of 
physics, e. g. the theory of elasticity and the theory of heat. 
For instance, if we combine case nr 2, which means a line with- 
out losses, with the boundary condition litt. B, we get a Green 
function of the form 


sl^X'shPls — 

Tjr , t \ U U 

K{x,^;p)=— — , x<i: 

P.sh'Ps 
u u 


K(x,g;p) 


sfp (s—x) sh^4 

u u ^ 

, 


u u 


.. ( 120 ) 
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which appears precisely in the study of the movement of a 
homogeneous elastic string with fixed ends or that of the air in 
an organ pipe, closed at both ends. The cases 2 and C combined 
give the Green function for a string with transversally freely 
movable ends or for a pipe with open ends. At last case nr 4, 
i. e. the cable line, where d=k sfp, has its counterpart in the 
theory of the flow of heat in a homogeneous slab of infinite 
length which is isolated to heat against the surrounding medium, 
its end being kept at a constant temperature. 

The expression (120) gives an example of a meromorphical 
basic integrator. In (119) again we meet a comparatively simple 
irrational integrator; this expression presents a many-valued 
function of p. 


Though the universal expressions (104) and (107) are very com- 
plicated, they can be governed comparatively easily by integratio- 
nal and operational methods if they be at first decomposed into 
simpler constituents. Such an analyzing procedure, the aim of 
which is to lead back the interpretation of a composed inte- 
grational expression to already known typical integrators, is 
characteristic of this branch of calculus. Usually, it will then 
be necessary to introduce expansions that are not closed, which 
means another kind of complication, but in many cases the state 
of things is such as to afford sufficiently good approximation 
by the use of only a few terms in the series occurring. 
Employing the coefficients of reflection 


Z+Z' 


and G" 


_Z^Z" 

'"Z+Z” 


( 121 ) 


we can easily transform the general expressions (104 a) for the 
voltage integrator into the form 


1 • p — 2 ^( 5 — ^ p—Sx 

-^rE',(p), 

l+'l 

z 

-F (®, p) = - 


l-G'G"-e-^^^ 1 + 


Z 


( 122 ) 
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In a similar way one can transform the expressions (107) for 
the Green function, getting 


K{xJ;p) 


e-S(i-x) 

26 r—aG''-e-^^^ 


when « < <f 


K{x,i;p) 


e-Xx-i) (i_G'.g-2«f) 
^2 6 


.. (123) 


when a: > * 


As we see, the same expression l—G'G" e~^^^ appears in all the 
denominators. In the most general case we have 


^^l/Vlp+l^i) (P+/'2). 

U' 


Now, in the complex integral (1) the path L between the points 
a ±ioo can be situated anywhere to the right of all singularities. 
We choose for d that branch of the square root which all over 
the path L has a positive real part. Then <; j along L. 

Further, if the reflection coefficients have suitable properties, 
e. g. if they are such that |6r'G"|<M when |p| is sufficiently 
large, one can remove L so far to the right that anj^where on 
this line 

|G'G"c-^^'«l<l (124) 

• 

Then the factor (1— can be expanded into a power 
series, and it is permitted to integrate term by term along L. 
By this means one gels the following series representations of 
the expressions (122) and (123): 

v\x, t) = W’ {x; p) • 6^':, (0, t) = W” {x; p) • A* (<), 


where 


W’(a;;p) = ~~ 
1 + 


Z' 

Z 


4 OO OO I 

lfc=u fc=«0 I 
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(125) 


W" (aj.p)= -^7- 

'+Z 


. [(2A;+1) s— a;]_j_ ^ Q'(k-hi ) . Q"k . J [(2A;+i) s+a:] (. 

A:=o k=o 

The series for the Green function becomes, if cc<^^: 

1 


K (x.^;p)= 




2d 


— S (2 Ars+I — x) 


k^-o 




^ — S i2ks~\~^'^x) 


I 


(126) 


A-=o 


- ^G’* • G''<* s-f-i] 

k~o 


<J(2Ar«— l+a?) I ^ 

A--=i I 


The expansion for the case tr>(f is obtained by making x and 
§ change places in this expression. 


The typical term in both the series (125) and (126) is of the 
form 

G{pHG'(p)V'-[G"{p)r e~^i (127) 

where ^ denotes a positive number and G (p) stands for one of 
the factors in front of the formulas. Accordingly, the solution 
of the original problem can be presented at once in an explicit 
form (we include, of course, in this term expressions involving 
quadratures with respect to Mn a finite number), if we already 
know the interpretation of the integrational expressions 

G{p)-H(t), G'(p)-H(t). G"(p)-H(i) and 
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Introducing the corresponding integral operators, we may write* 

cp{t)=G{pHG'{p)r-[G"{p)V'-e~^i-f(t) \ (j 28 ) 

=G X [G'] *=' X [G"] *” X f •/ (<) 


In the special case where the transmission line is distortion- 
less, i. e. d (p) = ^ (P+/0» we thus get, according to rule nr 7 of 
u 

our table, 


(Pit)=e •Gx[G']*'-'x[GT'X/ 


^--1 

uj 


(129) 


Here the symbolic operations ought to be performed in agreement 
with the general rules of Chapter 3. 


In the examples given above we have met some different types 
of integrational expressions. We will postpone the analogous 
treatment of the wave equation in three dimensions to a later 
chapter and now proceed to study more in detail some important 
classes of integrators in order to get general rules of interpreting 
the expressions already obtained and other in the practice fre- 
quently occurring types. 


7. A general theory of rational integrators. 

The expansion theorem. 

As already mentioned, integrators of this important class will 
always appear in connection with systems of ordinary linear 
differential equations with constant coefficients. Provided that 
the forces attacking have rational basic integrators, also the basic 
integrators of the unknown quantities will of course be rational 
functions of p, thus having no other singularities than a finite 
number of poles. 

Suppose F(p) to be such a function, regular at the point at 
infinity. In the neighbourhood of this point we have the ex- 
pansion 

CO 

f(p)=-2r 




* Cf. Appendix 6. 
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Further, suppose that F{p) has the poles po=0, pi, — pjsf of 
the orders fe, kv kn respectively. About the pole pn (n#0) 
an expansion of the following form exists: 

kn 

{F{P)} = -any+p-Pn ip-pn) (131) 

r‘=l 

(n = l,2***iV) 

Here Pn (p~^Pn) denotes a power series in (p—pn)- 

As to the expansion about the origin we prefer to write it in 
a somewhat different way, viz 

ko 

{FiP))o-^^p+p-Po{p) (132) 

r=o 

According to the definition of f (t)=F{p)-H (t) and to Jordan’s 
lemma, / (^) ee 0 for / < 0, and for ^ > 0 f{t) may be presented by 
a complex integral along a circle about the origin including all 
the poles of the integrand: 

= 

o 

t > 0 

Obviously, the explicit function / (/) will appear under different 
forms according to our starting from the expansion about the 
point at infinity or from those about the poles. 

P ( d) 

The residue of at the point oo is equal to minus the 

coefficient of p® in the series 
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Of course we must obtain the same result by applying term 
by term rule nr 1 of the table. 

Now, on the other hand, we will look at the poles inside the 
circle of integration. The residue of pn equals the coefficient of 

— - — in the expansion of 


P-Pn 


F (v)-—=e^n^-—‘F(p)‘e^^^n^K 

p p \rf 


Introducing the power series of — and we get for 

P 

the product of these functions : 

oo 


1 / t . , 

=pi.p,;-K+«--|i + • • • • 


Hence the residue of the point pn {n '4~ 0) l)ecomes 




For the pole at the origin we get the expansion: 


p» 

*/ A 


+P-Po(p)H 




and the corresponding residue is 
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Therefore, joining the contributions from all the poles we 
obtain the following form of the function f{t): 



ko N 

hn 

1 

/(<)= 

„,/=o n—\ 

^ 1 

Z"”” k-1 



•H{t) (135) 


In the practice one can usually arrive at this result in a 

F (v) 

simpler way, at first developing — - mio 'partial fractions 

P 

ing to well-known algebraical methods: 


F{p) 

P 


N 


kn 


LiP-pnY 


n—i v=i 


(136) 


Because F{p) is proposed to be regular at infinity we have 
no constant term in the first sum. Naturally we may as well 
directly make the disposition 


N 




'O^nv* 


(137) 


The coefficients a being determined, one has only to form the 
product F{p)’H{t) and then to apply term by term rules nrs 
1 and 3, thus obtaining the expression (135). 

This formula (135) in connection with (137) is often called 
the general expansion theorem. 

In the most common case F {p) has only simple poles, none 
of them at the origin. Then the expansion for F (p) becomes 


N 

J’(p)=ao+y^-an (138) 

L^P'—pn 

Tl=l 

and accordingly 



( 139 ) 
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Here the coefficients are very simply computed by the formulas : 

ao=F{0), an=—- lim {{p—pn)-F{p)} (140) 

Pn p=Pn 

If we write F (p) in the form of a quotient : 


J’(p) = 


’F(p)’ 


(141) 


and assume' ®(p) and ’F(p) to be polynomials in p, the order 
of ^^^(p) being at least as high as that of ®(p), the formula 
(139) will get the appearance: 


^p) 

^ip) 


H{t)^ 


N 


<p(0) I V <p( p«) 

Lpn-npn) 




(142) 


Hence, the explicit solution is immediately obtained, as soon 
as one has found the primary factors of the denominator. 

The expression just presented is of a simple form suitable for 
learning by heart; it is often named the Heaviside exyansion 
theorem. We will sometimes call it the simplified expansion 
theorem. 


In technical and physical applications F (p) will always be a 
function with real coefficients, and therefore all complex zeros 
of will appear in conjugate couples. U pn>,n"= — P ±i(o 

means such a pair of roots of ^l^{p)==0, we know that the corre- 
sponding terms in (142) must also be conjugate to each other. 
Hence, with the notation 

^ ^ 

their contribution will be 

{2A e-^' cos {(ot+a)}-n{t) (144) 

In all oscillatory mechanical or electrical systems not con- 
taining internal energy sources is positive on account of the 
losses. Thus, any couple of conjugate complex zeros of the de- 
nominator will give rise to a damped sinusoidal oscillation, setting 
in at the moment i^ = 0 with some initial value ( = 2^ cos a). 
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According to formula (134) the function in question is obtained 
in the form of a power series in t, which, as a rule, is Yery 
adequate for calculating / (t) for small values of t. On the other 
hand, the representation of /(t) by (135) has the advantage of 
giving a complete insight into the analytical nature of the function. 

Putting /=0 in (134) we get from (130) 

/ (0)=Ao=F (oo). 

Obviously, this important theorem, stating that the initial value 
of the time function is obtained, by putting p — oo m its basic 
integrator, always holds when F (p) is regular at infinity (thus 
even in cases where F (p) is not a rational function; we will 
come back to this question in later chapters).* 

More special, but often of practical importance, is the following 
theorem, which can be easily deduced from (135) and (137). 

If the rational integrator F (p) has all its poles situated in the 
region to the left of the imaginary axis of the p-plane, the sta- 
tionary value of the function f {t) — F{p)-H (t) is given by F (0). 


8. Extension of the expansion theorem to a class of 
meromorphical basic integrators. 

When F (p) is a meromorphical function, no singularities but 
poles will occur, and the expansion theorem can always be 
applied formally. Further, the correctness of such a procedure 
can be investigated by the well-known partial-fraction theorem 
due to Mittag-Leffler, In this chapter we shall make some 
examinations in that respect, confining ourselves for the sake 
of simplicity to the Heaviside form (139) of the expansion theorem. 

We will prescribe the function F (p) to have the following 
properties: 

1. Inside any finite region F (p) has only simple poles pi, p 2 * • *, 
which we assume to be numbered in such a way that 
I Pi I < I Pa I < I P3 1 < • The residues of F(p) corresponding 
to these poles may be ai, ch, as --. 

2. It shall be possible to determine an unlimited sequence of 
circles Ci, C 2 * • *, their centra being in {a, 0) where a denotes 


* Cf. Appendix 5. 
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a fixed and sufficiently large positive number, and the radius 
Bm tending to infinity with m, so that: 

a) none of them passes through a pole, 

b) F (p) is uniformly limited on all these circles, which means 
that |jF^(p)|< a for all points upon Cm, A being inde- 
pendent of m. 

3. The origin is a point of regularity for F (p). 

4. All poles of F (p) are situated to the left of a line L 
parallel to the imaginary axis of the p-plane. 


Now, P being any point inside Cm not coinciding with a pole 
Pr, according to the theorem of residues: 



the summation being extended over all the poles, in number 
n(m), enclosed by Cm. 

To begin with we transform the left hand side of (145) in the 
following way: 

_L flSv)..PL f 


f(o)+2| 


+Ar jip) .dp. 

^2niJ p ip-P) ^ 


r=l p 

Accordingly, (145) can be written 

n 

P 


F{P)=F{0)+S 




/^Py P—pr 2TciJ p {p—P) 

V=sl 


(146) 


The last term in this expression tends to zero when m (and 
thus also n) tends to infinity; indeed, its order of magnitude 
will be ^/Rm- Hence we have arrived at the expansion of F (p) 
in partial fractions (Mittag-Leffler’s theorem): 






(147) 
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The poles being all simple, the residue must be 

ay=lim [(p— (p)] (148) 

P^Py 

The theorem (147) is a direct generalization of (138). Putting 

^=ay we get even a formal agreement with the previous formulas 
Pv 

(138) and (140). 

It is easily verified that the rest-term in (146) tends uniform ly 
to zero when P being situated within a circle about the 

origin with a fixed but arbitrarily chosen radius q. This folio w 
from the inequality: 


Indeed, the right hand side will be < £ if m (or n) is made 
sufficiently large, and this independently of the position of P 
nside the circle Cq, 

In order to arrive from (147) at our expansion theorem we 
will now make use of the fourth condition prescribed for F (p). 
We denote by F a contour composed of a part of the line of 
iiitcgralion L in the de (initio n- formula (1) and a halFcirde, convex 
to the left, its centre being in (a, 0). Further we assume a number 
of N poles of F (p) to be situated kiside F, On account of the 
uniform convergence in (147) the integration may be performed 
term by term, and thus we obtain 



y^l 


Now, making the radius of the half-circle tend to infinity, we 
also cause N to do the same. Because of the second condition 
stated for F (p) it must be possible to draw the half-circle in 
such a way that | F (p) | <: A everywhere upon it; hence, by 


p{p-p) 


, Q‘ A' Pm 

{Ptn {Piii~ 


.. (149 
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Jordan’s lemma, its contribution to the contour integral must 
disappear in the limit, provided that < > 0, and therefore the 
integral along P has the same value as the integral along L in 
the definition formula. Accordingly, the expression (150) with 
jV=oo must represent the function F {p)-H (t) for positive values 
of t. Then, to be ready with the proof of the equality 




oo 







(151 


we only have to add that in the limit the contribution from a 
half-circle Cm, convex to the right, must be zero, if / < 0, and 
that the line L, on account of the 4th condition, may be deformed 
into such an infinite half-circle. 


We choose an example from Chapter 6 to elucidate the treatment 
of an integrator of the type in question. Thus, taking under 
consideration a homogeneous transmission line without losses, 
the ends being short-circuited, we will try to find the behaviour 
of the voltage along the line due to the initial state 

V (x), i{x,0)~0 (152) 

We assume no external fields to be present. 

According to (96), (98), (100), and (108) the expression of the 
voltage becomes 

v{x,t)^u-c-{r-\-pl)-J~G{x,^;p)-H{t)-j{s)-d^, ■■ (153) 


where the Green-function (Cf. (120)) has the following form 


shPx sJv {s — §) 

G (x, i;; p)= — — ® h 

shTs 

U 


sh^(s—x)sh^§ 


G (x, i;; p)= ; X > 

■sh"s 

u 


(154) 
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Thus our first task will be to examine an integrational ex- 
pression of the type 


®(P) 

W{p) 


Hit) 


sh^(s—x)-sh-i 

u 


shP s 

u 


.//(<). 


It is easily confirmed that this meromorphical integrator ful- 
fils the four conditions stated. The poles py are all situated 
upon the imaginary axis 

s 

Therefore we get, employing the expansion formula (142), 


^p) 

'F(p) 


■Hii) = 


0(0) 

>F(0) 


-f-OO 

, V’ 

y;sz OO 




(155) 


The apostrophe on the sign of summation denotes that the sub- 
script-value v—0 must be omitted. In our case the first term 
within the brackets disappears. Further, as 5 /i(m)=i’sin a and 
c/i(ia)=cos a, we get 

sinfi^*— (5— a:) j-sin 

I A A 

Py'''^'{Vy) i'Tt COS V7V (156) 

i . / ^ \ 

sin V* x\ 

V7t y s j 



Hence the expansion theorem gives the expression 


G(xJ;pyH 


oo 



y=l 



■Hit), 

(157) 


which is valid when x>§. 

In a quite similar way we get the function G {x, p)^H (t] 
when In this case, according to (154), we must obtain an 

expression identical with the one just found if x and § change 
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places; on the other side this expression (157) is symmetric in 
tT and^. Therefore we may conclude that (157) holds in any case. 

Now, inserting the results gained into (153) and then inter- 
preting the factor p as the symbol ", we get the following ex- 
plicit form of the voltage 


V {x, t)—U‘C» 



2 

7t 


8 CO 



0 v=l 




(158) 


The uniform convergence of the trigonometric series in the inter- 
val (0, s) permits us to integrate term by term. We denote by 
Uy the general Fourier coefficient of f(x), defined by 


ay 




(159) 


If j{x) be continuous, f {x) and f'{x) piecewise continuous in (0, 5), 
one can easily verify that v^-ay remains less than a fixed number 
when V tends to infinity. Therefore the new series converges 
uniformly with respect to t, and this will also be the case with 
the series obtained by differentiating it term by term. Hence 

it is permitted to put the symbol under the sign of summation, 

and thus we get the result: 



+ 


V {x, t) = 

oo 

y/ ? • * sin V sin ( y~ut 


CO 

\ Uy • sin I y—x 1 • cos ( y-ut 

^ \ ^ 1 \ « , 

j/=i 


■H(t) 


(160) 


Here and all the coefficients have the dimension of a 

voltage. 

E 



In passing we note that, the sense of the coefficients being another, 
the first sum in this expression represents the state of vibration 
of a freely movable elastic string with fixed ends starting from 
equilibrium position, the initial velocity being / (x). The last 
sum in (160), on the other side, gives the movement when / (it) 
denotes the initial position, the initial velocity being zero at all 
points of the string. Further we observe that, by a simple trans- 
formation of the trigonometric products, one gets v (x, t) presented 
as a superposition of an infinite number of sinusoidal waves 
travelling in both directions of the iT-axis with the velocity u. 


9. Irrational integrators in the form of powers of p. The 
notion of fractional integration. Treatment of certain inte- 
grators being rational functions of powers of p. 
Asymptotic expansions. 


In the second chapter we deduced the very important rule of 
calculus (18) 


P 




n 


being valid when n is a positive integer or zero. Now we are 
going to try to generalize this rule, investigating the sense of an 


integrator of the type — , s being any real number. For this 
purpose we will start from the integral of definition (1), giving 


(X-f loo 



O ’ — ioo 


(161) 


Because of the infinite extension of the line of integration we 
must set up the convergence condition 

5> — 1 (a) 

If 8 is not an integer, will be a many-valued function of p, 
the origin being a branch-point. We may, however, obtain a 
uniform branch of this function by making a cut from the origin 
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to infinity, conveniently along the negative real axis of the p- 
plane (Cf. fig. 4). 



The condition (a) being fulfilled, according to Jordan’s lemma 
we get the value zero, if ^<0, and if we take as path of in- 
tegration an infinite half-circle convex to the right. On the other 

9 

hand, if i > 0, we get the same result by integrating along an 
infinite half-circle convex to the left. 

Now, let L be an arbitrary loop starting and ending at the 
points — cxD ± 0 • i and enclosing the singularity at the origin. 
According to Cauchy’s theorem we may continuously deform the 
line L into the new path of integration (1 ->a) + + (^-►2) 
without influencing the value of the integral. The contributions 
from the quadrants (1— a) and (6—2), as already mentioned^ 
will disappear when / > 0. Hence, for positive values of we 
can use the loop L' instead of the line L on evaluating the 
integral. 

In order to lead back the integral along V to a wellknown 
function we will now deform the loop until it follows the borders 
of the cut closely, except at the origin, where it shall have the 
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form of a small circle. If, on this circle, we pul p=Q-e^ and 
\e^\ CO 1, we see that its contribution to the value of the integral 
will be of the order of magnitude q~~^. Hence, to gel a convergent 
integral when we must prescribe s to be negative. 

Accordingly, under the restrictions: 

- 1 < 5<0 {^) 


we may abstract from the small half-circle, thus getting the 
interpretation 


0 



pVt. 


dp+ 


27ti 




below the real axis above the real axis 


valid for positive values of L 

To examine this function of t and s we will occasionally put 

r = -^Sy 0<r< + l (}') 

Accordingly, 

0 CO 

— CO 0 

below the real axis above the real axis 


Now, in the cut p-plane we shall use that branch of the 
function 

pa _ • log p 

for which the logarithm has its principal value, i. e. 

7t > Im {log p) > — TT. 

Then, upon the lower border we get, putting \p\ = x, log p = 
log X — iTt, and upon the upper border log p = log x + ijt. The 
corresponding expressions for p^~^ are 

Therefore, as p=x-e^'^^=^—Xf we obtain 

0 CO 

■' ' ' 2m / 2m / 


0 
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or, combining the terms. 


oo oo 

/(<)=2^v J'x’^-^ er^ dx {\&2) 


This integral can readily be expressed by the so-called T- 
function of Legendre or the //-function of Gauss. Indeed, put' 
ting xt~u we get 



sin rTT 
rc 



Now, if Re{z} >0, we have by definition 


(163) 


oo 

Il(z—l)^r{z)= J'u^~^e~^du (164^) 

0 

Thus, instead of (162) we may write 

(165) 

By integration by parts in the definition formula above we 

find, provided that Ite { z) > 0, the recurrence rule 

/Z{s)=2-//{z-l) (166) 

Hence, n being a positive integer, 


n(z+n)={z+n)^lz+n-l)-^-- (z+l)-Jl(z) .... (167) 

Making si—O, as 11 (0) ~ F (1) = 1, we recover the well-known 
identity 

n(n)=n{n-l)--- 2 1==\n (168) 

heroin (166) or (167) we can obtain the analytic continuation 
of that function Tl {z) which by (164) has been defined in the 
region Re {z)> — ^. It will appear that U ( 2 ) is meromorphic in 
the whole ;s-plane and possesses the simple poles 2 = — 1, —2, 
— 3 - • Anywhere, save at these poles, the /I ( 2 )-function thus 
completed satisfies, besides (166), also the important relation 

ii(.).n(-.)=-^^ 

Sin TT 

as proved in textbooks on /'-functions. 


( 169 ) 
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Now, in the expression (165) just found for f (t) we will re- 
sin T 7t 1 

place the factor — — //(r— 1) by thus obtaining 

At last, returning to the notation s~—r, we arrive at the in- 
tegrational formula desired, viz 



The condition of convergence first stated, i. e. s> — 1, must 
still be sustained. The second condition, .9 < 0, on the contrary, 
may now be annihilated, because the critical integral around the 
small circle has disappeared from the expression. We understand 
that this condition must be immaterial already from the fact that 
it will not come in, if we integrate along a loop L' avoiding 
the origin; its introducing occasionally, however, made it easy 
to identify a constituent of the unknown function with the 11- 
function. 


Within the interval — 1<^<0 we have by (170) got defined 
a function of t not belonging to the general class of functions 
we have hitherto dealt with, as it tends to infinity when <-► +0, 

but obviously the expression — must be a basic integrator for 

'P 

this range of 5-values too. 

Because, if 5> — ■ 1, the function tends to zero 

uniformly with regard to arg p when p oo, we may for / > /o > 0 
employ the modified integral 

and here it is permitted to differentiate under the integral sign 
any number of times with respect to t. Hence, in the modified 
integral we have got a means of interpreting, for positive values 
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of U inlegrational expressions of the type 


p8—n 


H (t) (n=l, 2 ,--). 


But, s being any real number > — 1, we have obviously got a 

delinition of (t), valid for all real values of s. It is given 

by (170) if t>0; if ^<0, we postulate the corresponding time- 
function to vanish identically. It should be noticed that (170) 
changes into the more special formula (18), if n is a positive 
integer, and that the right hand side becomes zero for 5 = — 1 , 
— 2 •*, as the il-function gets infinite for these values of the 
argument. 

The integral operator corresponding to F(p)=p~~^ is 

t 

77(j;;.T A-')"-'* I (>’'> 


provided that 5 > 0 . Now, in Chapter 5 we have shown that an 
integer power of the symbol of integration 

t 

/=[ f it] 


should be formed according to the rule 

t 


I^ = 


1 


n—1 


{t — r) dr 


and it is natural to regard the previous formula (171) as a 
generalization of this result. Hence we are brought to the dejiniiion: 




f 

0 


(172) 


Such a power of the symbol I is sometimes called a fractional 
integration. Of course we need to prove that our definition is 
consistent with the addition theorem for common powers, viz 

(173) 

The proof is readily effected in the following way. 
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We form, starting from any function (p {t), the expression 

{rxp)x^{t)= 

t -1 p t 

^ Yl (s— ^9^(0 

0 J L_ 0 _ 

(174) 

t t 

0 T 

where the double-integral has been transformed according to 
Lejeune-Dirichlet’s method. Now, the integral 

t 

«/= J' ( 175 ) 

T 

is easily treated. Indeed, from simple substitutions results 
1 

“ (176) 

4 

on account of a fundamental theorem as to the il-function. 
Inserting (176) into (174), we find 'the following form for the 
symbolic product: 

t 



_ 0 _ 

But here, by definition, the right hand side is just 

By the same substitution as above, or directly, we get 

K 

n^fyy(<-0"-^-d^ xH(i)=-^^.H(<),(i78) 

- 0 

which formula must be identical with (170). 
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In the case where s<0 the symbol P will lead to divergent 
integrals. It is, however, possible to retain formally the definition 
(172) even in this case, provided that (p{t) is of the form 

(p(t) = 

In order to better elucidate this fact we will somewhat exa- 
mine the product PX(p (^), assuming 95 (t) developable into a 
Taylor-series 

00 

(ira) 

v=o 

Applying (172) to each term and making also use of the equality 

( 180 ) 

0 

as the right hand side has a meaning for all s and r, we may 
write for ^> 0 : 

00 

O"') 

If here we put 5—0, the left hand side will be iox(p(t)-(r{t), 
while the right hand side changes into the given series (179). 
If we assume 5 =n= a positive integer, (181) represents the result 
of an n-fold integration of ip (t) between the limits 0 and t. 
Further, putting s — —n, as il (— l) = i/ ( — 2 )= • • • = 00 , we have 

00 

r=n 

If, at last, we assume s — -~n—(jc, x being a positive number < 1, 
we obtain 


|/=0 

Here il {v-—n—x) will always be a finite number, and the expo* 
nent of t will be negative, as long as vcn+x. Thus, only 
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provided that all the derivatives of <p (t) of orders < «+« disappear 
at the moment t=0, the product I~^^X<p{t) can represent a 
function / (t) belonging to class nr 1 at least of the first chapter. 


In the foregoing we deduced the formula (170) by performing 
the complex integration in the definition formula (1). However, 
this rule may be more simply obtained by an opposite proceeding. 
In fact, computing according to (2) the function F{p), which corre- 


sponds to the time function 


■ 

n{sr 


we get 


F{p) 






(182) 


Here we put u=pT and assume, to begin with, that p is real 
and positive. Hence; 

CO 


1 




p* II (s) 


u^-er-^-du. 


Now, on account of (164), the last integral is identical with 
n (s), provided that s > — 1 . On the other side is an analytic 

r 

function of p, and this is also the case with the right hand side 
of (182) within a certain region. But, as these functions coincide 
upon the positive real axis, they must do so anywhere in their 
common domain of definition. Therefore they must also present 


the same analytic function, and we may accordingly put JP(p) = 


ps- 


The results obtained above furnish us with a means of inter- 
preting certain classes* of more complicated integrational ex- 
pressions, e. g. such as are rational functions of powers of p. 
Such an expression 


F{p) 



(183) 


* Cf. Appendix 5. 
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which we will first propose to be a basic integrator, Can be 
written more conveniently under the form 



1 +^ay-p—<’y 

1 -^Pr- 


(184) 


where s> — l and all Oy and gy>0. 

As before we make a cut in the p-plane along the negative 
real axis and replace the path of integration L by a loop enclosing 
the cut and wide enough that 



anywhere upon it. When this is the case, we may develop 
{1—U fiy-p~~Qyy~^ into a series and then, after multiplying by 
the numerator, perform the complex integral around the loop* 
term by term. Of course, we must arrive at the same result by 
putting at once 

<x> 

F iP)-H (0= " •(! + (t ) . . (185) 

.«=o 


and applying the general power-rule (170) to each term of this 
expression. The first term of / (t) becomes 

"■n(.)- 

All the others have an exponent >s. If s=0, we get /(0)=a, 
and it is easily seen that our previously stated rule (Cf. p. 60) 

F(oo)=/(0) 

subsists also for the class of integrators under consideration 
(Cf. Appendix 5). 

If in (184) — 1, jP (p) is not a basic integrator. The function 

P®-jP (p), however, satisfies the condition A of Chapter 2, and 
therefor^ we may employ the modified integral to get an inter- 
pretation of F {p)-H (f) for positive values of i also in this case. 

* Obviously the loop will enclos^/tjU^he poles occurring. 
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As a practical application we will examine the voltage transient 
on a coil with the resistance R and the inductance L, being suddenly 
switched to the end of an infinite cable in series with a battery 
with the resistance Ro and the electromotive force E. The impe- 


dance of such a line is Z' 


= i / and we get according to Ohm’s 
V pc 


law the following expression for the voltage: 

R"^pL 




Rq 


or 


v{t) = 


+B+pL+^‘ 


i+fi 

L P 


r 


E-E{t) (186) 


pc 


1 I 1 I ijr 1 1 

L p y c L pYp 




Hence the initial value is v (0)=E. Taking into account only the 
first four terms, we obtain: 

v{t) = 




V: 


r 1 
c L p^p 


1 




r 


rr 1 t% 
/cLliQ) 


, Ro {R-}-Ro) I 

’T7C I ‘ 


•]EH(t) 


(187) 


This expansion is fit for the calculation of v {t) for small values 
of L The general expression for v {t) is presented, as we know, 
by a rather complicated complex integral. The integrational method 
employed may be looked upon as a shorthand way of getting 
knowledge of the properties of this integral, when Ms a small 
number. 


If the origin is the only branch-point of the integrator F (p), 
the primary function / (t) contains, besides the contributions from 
the poles, if such are present, a constituent coming from the 
integration along the cut, viz 


i(t)= 


(»+) 




— CXD 


(188) 
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We will now briefly describe a method,* often applicable, of 
finding the analytical character of the function I (t) for great 
values of U 



If f is a great positive number, the exponential factor will 
obviously cause only the parts of the loop of integration nearest 
to the origin to give a considerable contribution to the value of the 
integral. We divide the loop into three parts with respect to a 
point — a on the negative real axis (Cf. fig. 5) and denote by 
Fu ip) the value of F (p) along the lower border of the cut and 
by Fo (p) its value along the upper border. On the former we 
have p=x er-^^=—x, on the latter = {x>0), and 

thus we may write: 

oo (0+) 

‘ W-i/ ^r)-FAv)]f‘d.+^^F^p)-‘-^-dp. . (189) 

a —a 


Further we assume that it is possible to find two real numbers 
K>0 and k, such that along the loop 

|jF(p)| < when x>a, (190) 

^ being any positive number. 

Then we can easily find an upper limit of the modulus of the 
first term in (189). In fact. 


A = 


\.F^iv)-F.m „ 


p~xt 


•dx 


(191) 


K , A K 7 IK ^ , 

< — / - •dx< / — - , if t>k. 

n J X Tt a f 7t a t — k 


* This method is further developed in Appendix 5. 
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Considering a, k, and K as given magnitudes, we may find by 
this formula a number h (s) such that, e being an arbitrarily small 
positive number, the inequality 

A<€ 

holds for all values of t greater than h. 

Now, to examine the behaviour of I (t) when t is great, we must 
further assume that in the neighbourhood of the origin F (p) 
can be written under the form 

n 

F {p)='^ay-p^y+p^ Rn{p) (192) 

V=1 


where 5 is a positive number and Si < <•• < Sn, and the 

rest-term such that, n being sufficiently large, 


\Rn{p)\<~e 

u 


(192 a) 


anywhere on the path of integration between the origin and the 
point —a. Then an upper limit of the contribution of the rest- 
term to / (f) is easily obtained. Indeed, 


(0 + ) 




s 

<- •£• 
a® 


(193) 


Here e may denote the same positive number as above. 

Thus, under the conditions stated, according to (189) and (193) 
we have, when t>ti: 


I 



< 2£, 


(194) 


To entertain the convergence when SyC 0, we may of course 
use a loop avoiding the origin. 

In the following we shall need an estimation of the magnitude 
of the difference 


( 0 +) ( 0 +) 

D= J^p^-e^^'dp— j^p^'C^^dp, 


( 195 ) 
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s being a real number. Putting, as above, on the upper 

border and p=x-e~^^ on the lower one, we get 


OD 

J) — — — 0—in8y 


from which follows 


oo 

\ D \<2 • J ' af - er -^- dx . 


(196) 


Now, if s>0, we have for all z > ^ the inequality a^< e*®. 


Hence, for such values of s , 

OO 


— , if ^ > • 

s e 


(197 a) 


III the case where ^<0, we put s=—r, thus obtaining 

oo oo 

\D\<2- /^~€r-^dx<~- fe~^'dx~^- \ ••• (197b) 

J J t 


We will now consider the difference 

( 0 +) 




y= 1 — oo 


( 5 +) 


v'' 1 r e**' 

y=l —a 

n (0+) (0+) 

^'p~ -e^^-dp 


v-l 



so 


According to (194): 



jy— 1 OO 



( 198 ) 


= 2 £+ 5 ^. 


Under the sign of summation expressions of the type (195) 
appear. Let us assume the exponents Si, 52, * * Sn to be such that 
for riL of them di> Su—1 > 0, 

for rh=n—ni of them ^^—l < —^2 < 0, 

and let a be an upper bound of the magnitudes of the coeffi- 
cients, i. e. \ay\<a, r=l,2*-* n. 

Then, in accordance with (197 a) and (197 b) we get the 
following estimation of S in (198): 


^ cc 

s <-■ 

7t 


o(<~) 


nv -- 
t- 




at 


di 

e 




(199) 


provided that ^ ^ 

By this inequality a number U («) can always be determined 
so that See whenever t> U- 

Thus we obtain from (198), taking t > the greatest of the 
numbers U and U just introduced (Gf. (191)), 


n ( 0 +) 

I dpi < 3 £ (200) 

P = 1 — OO 


€ being an arbitrarily small positive number. 

The integrals appearing in this inequality may be written 
symbolically in the form of integrational expressions, to be inter- 
preted when > 1 by the modified integral for positive values 
of t. Thus we have arrived at the required “asymptotic” form 
for the constituent I (t), valid when < is a great number, viz 


I {t)c^^^ay-p^yH (t). 


( 201 ) 
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According to (170) this means explicitly 

n 



y=l 

and we know that the error committed is less than 3 e, if t > ti, fc 
and the number of terms is sufficiently large (Cf. (192)). 


If the error for any fixed value of n tends to zero as t tends 
to infinity, we say that I {t) tends asymptotically to the right hand 
side of (202). It is easily proved that this may be the case if 
Si > 0. At first we remark that the contribution to I (t) from the 
rest term according to (193) is less than 


a"* 

fV-S 


1 (p) 1 


and therefore tends to zero with a. Now, let us assume a to 
approach zero as t increases, putting for example 


c 



On introducing this function into the inequality (199) we can 
infer at once that 0 as ^ oo. The same must be the case 
with the quantity A in formula (191), provided that the condition 
(190) be fulfilled in the whole interval x > 0, This, of course, 
necessitates that > 0. 

In order to get the function f {t)—F {p)'H {t) we must add to 
I {t) the contributions from occurring poles. If all the poles are 
situated in the left half-plane, their contributions must vanish 
with time, and because of (202) the stationary value of / {t) 
will be 

/(oo)=ai if 5i=0, /(oo)=0 if 5i > 0, and /(oo)=c)o if < 0. 

But in the neighbourhood of the origin the development (192) 
subsists, and thus we obtain 

f{oo)=F (0), 

which means an extension of the rule of page 60 to the class of 
irrational basic integrators in question (Cf. Appendix 5). 
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As a simple practical application of the method just described 
we will determine the asymptotic form of the voltage v (t) in the 
problem outlined above on page 76. For that aim we rewrite 
the expression (186) in the following more adequate way: 



1 +(B+J?o) • 




•E-ff (t) 

(203) 


Nearest about the origin we have the expansion 

F (p)=Ey^-{ie Yp-b (B+i?.)- ^^-p^- 
+ [L+B(B+B„)*.^]pl^+--j, 


and, as no poles are present, we get for great values of t 

v(t)oo (204) 






^ L+R{R+RoY ^ 


Vt 


2ts/t 


+ • 


We let aside, as a matter of pure calculus, the investigation of 
the domain of validity of this expression at a given degree of 
approximation, which is of course necessary in exact numerical 
calculation. 


10. Researches on certain special irrational integrators 
often occurring in practical applications. 

One of the first problems of great practical importance in 
electrotechnics concerned the transmission properties of a tele- 
graph cable. According to (116) and (119) the behaviour of the 
voltage along a cable of infinite length is determined by the 
equation 

zr— v/rc-v/p*® 

— e{t) 

l+S/l-Vp-Z’ip) 


V {x, t) 


(205) 
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In many cases the impedance Z\p) of the transmitter can be 
regarded as constant with respect to p, and often its value is so 
small that one can equal to 1 the denominator of (205). 

Thus we are confronted in this problem with an irrational 
integrator of the type a being a positive number, and our 

nearest task will be to find the corresponding primary function 


<p{t)=0{p)-H{t)=e-“^-H{t) (206) 

According to the integral of definition (1) 

G+ioo 



<r — too 


To get a one- valued integrand we make a cut along the negative 
real axis (Cf. fig. 4). Putting where — 

we have 

er-^\/p==e^\/Q-{ooB~^+i • 8in|) . 

From the equality | ® (p) | 2 it follows that 

tends to zero as | p | -► «>, uniformly in the angular region stated. 
By Jordan’s lemma we conclude in a wellknown way that9)(t) = 0 
when t< 0, und that for t>0 the line of integration L may be 
replaced by a loop L' enclosing the cut. Therefore 

t>0 (208) 

v 

Differentiating under the sign of integration, we obtain an integral 
that converges uniformly along the loop, if t > to > 0. The inte- 
grand being finite at the origin, we may contract the loop until 
it encloses tightly the cut and then present 9>' (f) by the integral 

( 0 +) 

(P r^\/i+Pf-dp, < > to > 0 (209) 

— 00 

According to the identity 

pt-«l/p=t.[(l/^-^)-^,] (210) 


<P(P) 

P 
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this can be rewritten 


( 0 +) 

.J_. 

27tiJ 




■dp. 


< > <0 > 0 . 

Now we will employ another variable y, putting 
Vtp—^=iy, hence p=\‘iiy 


2)/t 


2 )/< 


and 


dp-=—-liy-{ “ 


2^1 


dy 


(2111 


( 212 ) 


As to the image of the path of integration in the complex ?/- 
plane we note that: 

When p moves from — oo to 0 , t/ moves from — 004. to 

2|/t 2}/t 

When p moves from 0 to —00 y moves from to 00 + 

2j/t 2Vt 

Consequently the path of integration in the j/-plane must be a 
line yl, parallel to the real axis (Cf. tig 6). 



The transformed integral looks: 

A 


f > > 0 


Here the integrand is holomorphical within any region bounded 
by the real axis, the line A, and the segments between them 
of two lines, parallel to the imaginary axis; as obviously the 
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contributions from the latter boundaries disappear, when their 
distances from the origin tend to infinity, we conclude that we 
may take the real axis as path of integration instead of A. 
Further it is readily seen that the first term under the sign of 
integration in (213), viz i y-e~y\ being an odd function, will 
give no contribution to the value of the integral. Hence the 
expression is simplified to 


fit) 



■ CO 


... (214) 


t>to:>0 


We are now in a position to introduce the so-called error' 
functioriy which we will present in the following form, employed 
by many authors: 


Obviously 


hence 


X 



Here the right hand side equals cp' {t) but for the sign, 
we have 


Thus 

(216) 


Integrating on both sides between the limits 0 and t, we get, 
as erf {oo) = i, 

(<)==! -erf |-~J, f >0 (217) 


Accordingly we have won the new inlegrational rule: 


( 218 ) 
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By (214), as 9 j( 0)=0, the corresponding integral operator becomes 


<D= 



o* 

4«— t). 


dr 

{t-ryi* 


(219) 


For the rest, formula (214) provides us with another important 

( 220 ) 

2)/7r t^t 

Thus we see that also the expression p has the character 

of a basic integrator, and the corresponding primary function 
vanishes for t—0. 

Differentiating both members of (218) with respect to a, we 
obtain one more practical rule of calculus: 

(<) » H (t) (221) 


The integrator does not admit an expansion into a series 

of powers of However, one can easily derive an asymptotic 

expression for (p (t), valid for small values of t. Indeed, by 
repeated integration by parts in the definition formula (215) one 
finds, as erf (oo) = l, 

1— erf (*) = 


3 , 

2 2 


• ar“'— • - ■ 

f/yr L 2 2 


( 222 ) 


the rest-term being; 


1_ 1 3 (2n— 1) 

]/^ 2 2 




2n+i 


■dp 


n > 2 


a<> 


(222 a) 


It is readily proved that 


|lJ2»--xl< f •••• (222 b) 

fr TT 2 2 2 



87 


From this follows that the rest-term is always numerically less 
than the last term of the precedent expansion. Thus we have 
got an asymptotic representation of 1 — erf {x), available if a? is a 
great number. 

If specially we insert obtain the wanted asymp- 

totic expansion for cp (/), viz 

5l)(I) = l- (223) 

Vrc a L 2 2 2 \o?! 2 2 2 Wl J 

suitable for small values of t or, more strictly speaking, of ^1, 

X 

as we have a—X‘^rc, and the error being less than the last term 
taken into account. 


Another asymptotic expression for the function 
cp {t). 


applicable for great values of t, can be deduced by the general 
method in the foregoing chapter. In fact, upon a loop close to 
the cut (Cf. fig. 5) = and thus the contribution to the 

integral from the part of {he loop to the left of a point —a 
will be arbitrarily small if t is great enough. On the other side, 
this point can be chosen as near to the origin that the function 
0 -a\ip jjg approximated by a number of terms in the begin- 
ning of the corresponding power-series in aj/p. At last, the path 
of integration can again be completed by the parts to the left of 
the point —a. Obviously, the result must be the same if we 
expand formally the integrational expression, putting 


r N 


(pit) 


I 




Iv 


■Hit). 


(224) 


and then, when t> 0 , interpret each term hy the modified inte- 
gral. Hence we get, according to (170), 
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9 


N 






(225) 


Here all terms corresponding to even values of r (v=0 excepted) 
disappear, and the first part of the expansion runs 


(P(t) <^3 1- 


T-fL 




y-n:Y2^i 3 \2|/</ |2 5\2V'f 




(225 a) 


Referring to the general investigations in Chapter 9, we pass 
over a more detailed examination of the procedure. 


In Chapter 6 has been shown that the treatment of the uni- 
versal equation of telegraphy leads to integrational expressions of 

the form where (Cf. (126)) d~- (p+f-ii) This 

u U 

expression can readily be led back to an integrator of a more 
symmetric structure, the primary function of which we are now 
going to determine. Hence we introduce 

Z(f) = -7=£= e-“\/^^^.ff(i) (226) 

assuming a positive and employing that branch of the root which 
tends to p when p tends to infinityi 

If we rewrite the expression in the following way: 


X (t-a) (227) 

YP^ + O)^ 

we infer that x(t)=0 when t<a, and further obtain a new 
integrator 

® iv ) = (228) 


which satisfies the conditions for Jordan’s lemma and obviously 
is a basic integrator. 

We will deduce the explicit time function in two different ways. 
At first we develop the exponential function into a series and 
treat term by term the integrational expression thus appearing. 



This is permitted, since ©(p) is regular at infinity, as we have 
chosen the branch of the root. Accordingly we get as a first 
result: 

x(0= 


I 


( — ao)Y p — pV 




0 )^ 


0) 


H (t-a) . . (229) 


v=0 


Here we are faced with another important class of basic inte 
grators, which may be defined, more generally, by 


U{p) 


JP 


I ^ /’ 


(230> 


the exponent / now being any real number > 0. When i is 
not an integer, there exists also a branch-point at infinity. To 
get a one-valued function we make a cut consisting of two parts, 
one between the branch-points ± i w and another along the 
negative real axis, as shown by fig. 7. For ^p^+co^ we choose 
as above that branch which is positive when p is real, and for 
()/p^+w^—p)^= exp U*log p)] the branch correspon- 

ding to the so-called principal value of the logarithm, i. e. that 
which remains real when p is real. 



Fig. 7. 


When p -> oo in the cut plane, JJ (p) tends to zero, if > 0,. 
but to 1, if x=0. Obviously 17 (p) is a basic integrator; the 
primary function u {t) ^ 0, when f < 0. When t > 0, employing. 
Jordan’s lemma, we may define u (t) by a complex integral 
extended along a loop L\ tightly enclosing the cut (Cf. fig. 7). 
By the change of variable of integration p — ojz, which geometri- 
cally only means a uniform reproduction of the path, we get 
the expression 
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u(t) = V{p)-H{t) 


= J_ /'(I' 
2iiiJ 

V 


>+l-z)^ 


j/z^+l 


•e“"-*-dz 


(231) 


t >0 


Here the right hand side just equals a well-known integral 
definition (which holds even when X is complex) of the so-called 
Bessel function of the order X and the argument (x)t. Thus, em- 
ploying a common notation, we have 

u{t)=Jx{o}t), when <>0 (232) 

In the special case under consideration (Cf. (229)) only integer 
values of X occur. Then there is not a singular point at infinity, 
and the cut along the negative real axis can be omitted. 

Applying our new integrational rule 

\X 

^Vh {t)=Ji i(ot)-H {t); A> 0 (233) 

to the expansion (229), we get the result 


p lyp^o)^-— 

I (0 




fp^+ 




or 




{—a(oy 


(234) 


We won’t discuss more in detail this form of the solution. By 
introducing into (234) the power-series of Jv {(o{t—a)) (>^=0, 1, 2* • •), 
taken from textbooks of Bessel functions, one shall get a power- 
series in o)(t—-a) for % (<), suitable for numerical evaluation of the 
function, when oj(i—a) is small. Of course, one can also obtain 

the same expansion for x (0 ^7 developing 0 (p) into powers of - 

'P 

and then applying the first rule of the collection in Chapter 5. 


The second mode of obtaining the function x {t) is based upon 
certain transformations of the definition integral 


X 


<r+ioo 





(235) 
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If t> a, we may take as path of integration a loop enclosing 
the cut between the points ± io) (Cf. fig. 7). We also prefer 
another variable of integration <p, according to the equation 

p=iwsin9> (236) 

When p moves along the right border of the cut from — i o) 

7t TC 

to +i 0 ), (p will pass from — — to + „ through real values, and 

^ A 

again, when p moves along the left border from io} io —i o), (p 

77/ 7t 

will return from — to — — . Further 

A A 


±|/"pM^^=co cos (p. 


On the right side of the cut arg (p^+w^) =0, thus /p^+w'^ > 0; 

on the left side arg (p^+w^) = 27r and >/p^+^ < 0. At last we 
d/D 

have i‘d(p, the + sign corresponding to the right 

yv^+(o^ 


border, the — sign to the left one. Therefore we get instead 
of (235), 


X(t) = / • sinf-aoi • cos ff . 

27r ' 


27t 


' ^iw^singp'+awcosy' . 


TT 

2 


TT 

2 


or, putting (p—7t—(p in the last integral, 


~2 


^ (/) = - 1 — j * sin^r— ottj • 00S9P . _J_ j ^iuit • sm(f— aio • coBf . 

27^ 27r / 


IT 

2 


TT 

2 


Here the integrals may be joint to one, and as an integration 
between the limits — ^ and ^ because of the periodicity equals 

A A 

an integration between 0 and 27r, we obtain for t > a 


X{t) 


27 T 

-2nJ 


sinT’—ooi • oos9t> . 


(237) 
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Into this expression we will introduce an auxiliary imaginary 
angle 9, defined by the equations 

— a^'sin 0, — a^ cos 0. 

Then icot- sin (p—ao) - cos (p=io)y^t^—a^^cos {(p—9) and accordingly 

f (238) 


27V J 



But here we may replace rp — 0 by (p, retaining the limits and 
the path of integration. Indeed, if we put q) — 9==(p\ the new 

path between 9)'== — 0 and (p'==27t—9 
should be represented by the straight 
line ^ B in fig. 8. Because of the 
periodicity of the cos-function an 
integration along 0 A gives the same 
result as an integration along 27r — B. 
Further, on account of the regula- 
rity of the integrand we must get 
the value zero on integrating around 
the parallelogram 0—^— JB— 27r 

— 0. Therefore, the result must be 
the same if one integrates along A B or along 0 — 27r. 

Thus we have obtained the required function x (t) in the form 
of a definite integral: 

27r 4 

(<-«) (239) 

0 

Now the Bessel function of order zero can also be presented 
under a similar form, viz 

*2n. 

(240) 

27tJ 


and thus our result can be expressed by the equality 
X (<)= j, (i-a). 


(241) 


+{/)•“ 

In the following chapter we shall get use of this formula when 
studying the general solution of the equation of telegraphy. 
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As a last example of methods of treating irrational integrators 

V 

we choose some researches on the expression which 

often appears e. g. in the theory of electric filters and recurrent 
networks. As this corresponds to the special case i=0 in (233), 
we know d priori that the function 

(242) 


when t > 0, coincides with the Bessel function Jo {(lit). 

This basic integrator, because of its simplicity, is fit for illu- 
strating the above mentioned method of developing F (p) into 


powers of 


Thus, writing 





and expanding in accordance with the binomial formula: 


-1 


2 

r 


L< W' 

i6 jxy 


\2/ (1 12/ 

(|3)^‘\2/ 

'■(O? 



(243) 


we get on employing term by term the rule 
the result: 


• «(<)= 


1 

■H {t)=^-H{t) 

p” ' ' ' ' 


(244) 


[ 






Here the brackets just enclose the well-known series for Jo (cot). 

The correctness of the operations performed is easily verified, 
if in the definition integral for > 0 we integrate e. g. along a 
circle about the origin, its radius being > co. 


Obviously, the function u{t) to be determined will appear in 
varying analytical forms according to the path of integration 
chosen. If we integrate along the loop enclosing the cut between 
the points ± iio (Cf. fig. 7), we get by employing the variable 
V 

x=-~-, when t > 0, 
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l 

2 /" 1 

u{t)==~y ( 246 ) 

0 

This definite integral, as well as the one obtained by putting 
tr=cos z, i. e. 

7t 

a 

u{t) = -- /cos [{wt)- cos z]-dz (245 a) 

0 

represents a frequently used definition formula for Jo {cot). (Cf. (240)). 


Now we are going to examine more in detail another integral 
form for u {t), which we obtain by using two cuts between ± i(jt> 
and oo, parallel to the real axis. (Cf. fig. 9.) 



For positive values of t we may integrate along the loops 
L' and L'\ tightly enclosing these cuts; next to the branch-points 
we first assume the loops to have the form of small circles; 
however, it is easily proved that the contributions from these 
circles vanish when the radii tend to zero, though the integrand 
becomes infinite at the points ± uo. 

Upon the lower loop we put p=— — iw. Hence 

^ % • arc tg 

Upon the upper loop we put p=— Hence 

. , — 2(iJ 

, farotg -g- 

i-2/Sw=/?-|//e»+4w»-e ^ . 
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Occasionally we will use the abbreviations 

- {(}), - ^-arc tg (/g) (246) 

Then the expression under the integral sign in the definition 
formula 


1 


ppt 


2mJ 

VL" 


(247) 


gets the following form on the different parts of the loop: 
(— oo I — ico) -► ( 0 I — iw) : +-4 

( 0 I — jw) -► (— oo I — ia>) : — ^ 

(— oo| 0 I i<D) : 

( 0 I iw)->(— oo| io}) : 


Combining the different contributions, we therefore obtain: 



(yS) sin [wf-a {§)ld^ .... (248) 


Hence, our original complex integral (247) now has been 
transformed into a real integral, and we get for t > 0 after sub- 
stituting the expressions (246): 


•g~^ sin (wf-f^arc tg (249) 


According to well-known trigonometric formulas this expression 
can also be written 

u{t) = 


}/2 / V'/8+V'/J"+4w^ 


e~^^-dp-sin od 


^ J 1//9Y/5«-|-4w" 

0 

C 50 

/ 7 - ^ — ■ ■ • e~^ • dp • cos lot. 

^ J )//J-|//5“-H4w“ 


. . . ( 250 ) 


0 
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From this representation we can easily deduce an asymptotic 
expression for u{t), available for large values of t If Ms a 
great positive number, the factor even for small values of 
will be so minute that the integrand is practically zero. Next to 

the origin the factor is dominating. On the other side, the 



factor ^±^+yp^+4tn)^ is continuous at the origin, 

where it equals y=' reason we might replace the given 

integrand by the simpler one 

1 

)/^V^ 


retaining the limits 0 and oo. Then we shall obtain a value of the 
right hand side that must agree better with the true value the 
larger t is, thus furnishing us with an asymptotic expression for 
u{t), viz 

u{t)oo — {t)-{sin cot+cos cot), (251) 

TT- j/w 


where 


oo 



0 


(251 a) 


By the substitution ^ the, last integral is essentially 


identified with the so-called error’integral, and we find 



Therefore : 



(252) 


Here we recognize on the right hand side a well-known 
asymptotic form for the function Jo (cot) (Cf. e, g. Funktionen- 
tafeln, von Jahnke-Emde.) 

Of course, by developing the irrational factor under the sign 
of integration into a power series of one can arrive at more 
complete asymptotic expressions for u (t). This shall be done in 
a somewhat different way in Section 5 of the Appendix. 
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11. Applications to selected problems of physics and 
electrotechnics. 

In this chapter we are going to employ integrational and ope- 
rational methods for treating some important problems of theo- 
retical physics and electrotechnics, and we will generally try to 
work out the solutions in such forms as are already more or less 
familiar, the comparison of the new methods with older ones 
thus being facilitated. 


A. The differential equation of a freely movable homogeneous 

elastic string. 


If u denotes the elongation, a a constant of the dimension of 
a velocity, this equation is 


^ dx^ 


(253) 


We shall now integrate the equation under the initial conditions 

fini 

u—F{x), —=G{x), when <=0 (254) 

oz 


For that aim we introduce a basic integrator depending on x, 
writing 

u{x,t) = U {x; p)-H {t) (255) 


Then, according to the differentiation rule: 
•!^=[p.U-p-F{x)]-jH{t) 

U-jf^-F {x)~p-G {x)]-H it) 


(256) 


(Cf. the Note on page 45). 


Inserting (256) in (253), we obtain 

[p 2 . U-p^-F {x)-p-G {x)]-H (t)=a“ 


d^U 

da? 


■Hit). 


Hence, according to the rule of uniqueness we get for the basic 
integrator the ordinary differential equation: 


d^U 

dx^ 






(; 


(257) 
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The corresponding homogeneous equation having the two in- 
dependent solutions 

and e “ , 

we may put, employing the method of varying the constants, 


(258) 

With the occasional notation 

k{x) = -^-F{x)-^-G{ic) (259) 


we get 


A'{x)=l-e “*•&(«), = .... (260) 

^ A 


Hence the required integral of (257) will be 
1 ^x r-^i 1 -^x 

I7(®;p)=^-e“ • / e “•.A;(.^).d^-^e ^ ■ e “ -ki^ydi. (261) 


Making both the members operate on H (t), we find 


«(«, t) =— - /e “ •H{t)-k{i)-di—^ J e 


Hit)-k{i)-di 


But, as we know 




and further 


= 0 if x+at; H\t 
Hence we may write as well 


_ n . 


= 0 if ^ < x—at. 


u{x,t)=-\- fk{g)-d§-\- fk{§yd§. 
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Inserting k (^) from (259) and combining the integrals, we get 

x+cU X'\-(xi 

u {X, Jf f ? .... (263) 

Xr—oi X — di 

Here p is equivalent to the symbol , and we obtain according 
to the identity 

</ia) 

Tf (x)-dx=q>’ {a)-F {q> (a))-/' {a)-F (/ (a)) (264) 

f(a) 

the final expression for the elongation of the string: 

x+at 

« (®> 0=^ [F (ic+a^)+i^ S {s)-ds • • (265) 

X — at 

This is easily proved to be the only solution of (253) satis- 
fying the initial conditions prescribed. The expression (265) 
represents the so-called d'Alembert's form of the solution. 


B. The general wave equation in three dimensions. 

The homogeneous equation looks: 

— =V.A«, (266) 

a being as before a velocity, and A denoting the Laplacian 
symbol. We will assume the initial conditions 

(«)<=»=F (a;, j/, z), (267) 

Putting u (*r, y, z,t) = U {x, y, z; p)'H (i), we arrive as above at the 
equality 

[p2- F {x, y, z)—p^G {x, t/, z)\H (t)=a^‘ A U H (t), 

which proves that the basic integrator U {x, i/, z; p) must satisfy 
the partial differential equation 

A U-'^~V^-^,-F(x,y,z)-^^.G{x,y,z). 


(268) 
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Now, it is a well-known fact, that the equation 

AV-^-V=-4?r-9(®,t/,z) (269) 


has particular integrals of the form 

y(®,y,z)= (270) 

doj being an element of volume, and the integration being extended 
over the whole space. (If the Oc^nsity of mass q extends to in- 
finity, it must be prescribed to decrease sufficiently rapidly with 
an increasing radius-vector r). 

Hence, in our case (268): 

Q{x,y,z)=^--'^^-F{x,y,z)+-j-^^-G{x,y,z). .. (271) 


For interpreting a volume integral of the type 

i) 


I {X, y, z, t)^ ff f^ ° H {t)-dw .... (272) 


we conveniently employ polar coordinates (r, 0, 99). The pole 
being situated at (a?, y, z) we have the relations : 


Accordingly, 


^=a7H-r sin 0 cos 99, 
7 ]=y+r sm 0 sin 99, 
^=2;+r cos 0, 
dio=r^'Sin d dr-d9 d(p.} 

v{x, y, z, t) = 


(273) 


00 TT 2 Tt 

/// (iT+r sin 0 COS 99 , y+r sin 0 sin 99 , 2+r-cos 0)* 

000 


•IT ±^j-r*sin 0-dr-d0^d9P (274) 


The factor H differs from zero and equals 1 only when 

d: r > 0, i. e. when T r < at. Now, at and r being positive, the 
condition — r < at is eo ipso satisfied. The second condition. 
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r < at, brings a limitation as to the variation of r. Obviously, 
instead of integrating with respect to r between the limits 0 and 
oo we may as well integrate from 0 to at, thus obtaining 


v{x,y,z,t) = 

at IT 'In 

•Q {x+r-sin d coscp, i/+r sin 0 sin <jP, z+r cos 6)- 

sin 9-d9-d(p (275) 

In the expression (271) also a factor p (or p^) equivalent to 
a differentiation symbol occurs. By (275) 



p-v{x,y,z,t)= (276) 

TT 271 

a?-t- J' J'^ (a:+ai sin 0-cos (p, t/H-ot-sin Q -sin 90, z+ot-cos 0)- 

•sin 0-d0-d9). 


Hence our solution of (266) will appear in the form 


d 

dt\ 


n 271 

Wf’^ 


u{x, y, z, t) - 


(277) 


(a; + at- sin 0- cos 93, t/ + at- sin 0- sin 90, z-Hot-cos 0)- 


•sin 0-d0 


n 2n 


■#] 


+• 


or 




G {x+at sin d'co^ q), y+at sin d sin cp, z-\-at cos 0)* 
•sin O-Ud’dfp, 


u{x,y,z,t) = -[t-F{at)]+t-G{at), (277 a) 


if we denote by F (at) the mean value of the function F (x^ y, z) 
upon a sphere with radius at about {x, y, z). 

Thus we have arrived at a well-known form of the solution of 
our problem and may refer to textbooks of mechanics as to the 
discussion of its physical significance. 

In the case of an inhomogeneous wave-equation, which may 
be written 

I d'^u 

= A u+iTCQ {x, y, z, t). 


(266 a) 
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we get instead of (268) the integrational equation 

A17 {X, y, z)-^^-G(x, y, z)-4?r -B (®, y, z-, p), (268 a) 

Of d Of 


B {x, y, z; p) denoting the basic integrator of the density function 
Q (x, y, z, t). Accordingly, to the solution of the homogeneous 
equation just obtained we now have to add the term 


Ml {x, y, z, t) 


-Iff- 




(278) 


But here 

B (I. p)‘e^~a H (<)=e*o^^ {§, rj, ri, i;, t ±^j- 

Now we will make the assumption that also the volume den- 
sity, being a function of time, keeps a constant value until the 
moment <=0. In that case the solution corresponding to the -|- 
sign obviously must be omitted, and we get 


ui{x,y, z, t)= 


m- 




T ' 

— l*dw. 
a 


(279) 


Here the integral expresses the so-called retarded 'potential at 
the point {x^ i/, z) under consideration. 


C. The equation of conduction of heat in three dimensions. 

If we denote by a a constant of the dimension L*T—^ the 
homogeneous differential equation in question has the appearance 


(280) 

We shall try to find a solution that for t=0 coincides with a 
given function of space, say 

(u)t^o=F (x, t/y z) (281) 

Introducing u (x, y, Zyt) = U {x, y, z; p)'H (t), one gets the follow- 
ing equation for the basic integrator 

AU-l,-U=-P--F{x,y,z) (282) 

u a 



According to (269) and (270) this equation has particular so- 
lutions of the form 

V {X, y, z; p) . . (283) 

We will employ the notation: 


(p{r, t)=e^a -Hit).. 


thus obtaining, when t> 0, 


u (x, y, z, t) = 4 ^ f 

The basic integrator appearing in (284) is of a type already 
treated in the foregoing chapter. Thus, by formula (2l8) 


for the applicable branch of the square root, and accordingly, 
when t> 0: 



Therefore, in agreement with (285), the solution of our problem 
must be 

or, if we again employ polar coordinates (r, 0, y) and, to sim- 
plify our mode of writing, introduce the mean value of F 

n Stt 

F{r)=^ffF{i,v,S)-^ine-ded<p (289) 

0 0 

upon a sphere of radius r about (xyz): 


u (*. y . 0 - f\r)-e ^-r^ dr. 


But this is just the classical diffusion-formula. 
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If we write the inhomogeneous equation of conduction of heat 

Au+in:Q{x,y,z,t) (280a) 

we get in this case, instead of (282), the integralional equation 

U = y, z)-^ R {x, y, z; p) (282 a) 

CL d CL 

and accordinglj^ we have to add to the above solution (283) a 
solution TJi of the equation 

Ui= R {x, y, z; p). 

Cl Cl 


This solution is obtained from (270), giving 




Attaching the unit-function H (t) to both members, we get on 
account of (287), as ^(r,0)=0, the following expression for the 
complementary function 

Ui (x, y, z, 

t 


Here, if suitable, the integrations Vith respect to w and r may 
generally be interchanged. 


D. The universal homogeneous equation of telegraphy. 

In the universal form of the wave equation time derivatives 
of the first as well as of the second order appear. In the one- 
dimensional case the equation of telegraphy offers a practically 
important example of such an equation. According to Chapter 6, 
formula (94), where we shall put B and E equal to zero, the 
voltage as well as the current is governed by the equation: 




(292) 
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We will now integrate this equation under the initial con- 
ditions 

V {X, 0)=F (x), ^{x, 0)=G {x) (293) 

01 / 

If we put V (<r, t) — V{x; p)'H (t) and introduce 6 (p), defined 
by (99), we obtain: 

^-d\v)-V=-p[vdc-F{x)+{al+rcyF{x)+lc-G{x)l{2U) 

Note. From an electrotechnical point of view it is more natural to assume 
the initial conditions 

v{x,0)=F{x), i{x,0)=K (x) (293a) 

Then the integrator equation looks 

dW 

= [c{r+pl)-F(x)-l-K'(x)] .... (294a) 

This problem is therefore analytically equivalent to the former. 

Accordingly, we have to treat an equation of Ihe type: 
dW 

^^--d%p)-V=p-f{x) (295) 

Making the attempt 

V(x; p)—A {x)-e^^+B {x)-e~^^, 
we get, employing the method of varying the constants, 

-4'(a:)=— • e-^^-j{x), B' = {x) 

and from this: 


F(®;p)=^- ye-‘5(--f)-/(i?)-d/. .. (296) 

a 

Accordingly, the problem under consideration is essentially 
solved, if we know how to interpret the special integrational 

expression (.C>0)* We have already encountered this 

expression in Chapter 6, and now we will determine the corre- 
sponding primary function 

( 297 ) 
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For that scope we conveniently rewrite the transmission expo 
nent d{p) under the form 




u 


where 


1 Ir , a 


1 Ir a 


\lc 


2\l c 


1 


" 2 ' c 

Hence we can write x (^. t) more completely 




V (p+kY — 


e « 


Hit). 


(298) 

(299) 

(300) 


We may easily deduce the explicit function by making use 
of the following integrational rule, which has been proved in 
the foregoing chapter: 

P — • (f)=j, (t-a). .. (301) 

y p^+o>^ 

If we put o)=iA and employ the common notation J<,{ix)= 
lo (x), we get another appearance of this rule, viz 


P 


\lp^-?? 


jj (t^a) . . (301 a> 


To make this formula applicable to our problem we must at 
first somewhat transform (300) acgording to the so-called dis- 
placement rule 

O ip+K)-f {t)='e~’'‘.0 (p)-e’’* f (t). 

Thus we get, writing 


X(S>t)-U‘ 


P + K 


y/ {p+Kf—A^ 


—•e w 


- JL 

P + K 


Hit), 


because 




’’■Hit) (302) 

Vp*— 

Now (301a) can be directly applied, and we obtain xis,t) 
der the. form 


un- 
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( 303 ) 

After that we will return to formula (296) and make both 
members operate on H (t), thus getting 

a X 

v{x,t)=-\J x(§-x,t)i{§)-d§-\Jx{x-§,t)-i{§)-d§. (304) 

X /? 

Now, by (303), x{§—x,i)^Q if §>x-\-ut, and accordingly the 
limits of the first integral may be as well x and x-\-ut. Further 
X{x — t) = 0 if i<x—ut, and thus we can take the limits x — ut 
and X in the last integral. When the function is not zero, we 
have also Z (— ^, t)=Z(i, <). Therefore the last equality (304) 
may be written 

x+tU 

v{x,t)=—^^- J' K{x-§,t)-i{i)-d§ (305) 

From this we obtain by well-known methods 


—^•[x{—ut,t)-f(x+tU)+ 7 ,{ut,t)-f{x—ut)\ (306) 

x+ut 

-l-J'fi^)-j^X{x-i.t)d§. 

X — tU 

Further it follows from (303), when ^=— «f+0 or 0, that 

Xi±ut,t)=u-e^‘‘*Io{0)-H{arg>+0)=u-e^’^, (307) 

and from (305) and (306) we get 

x^vJt 

r; {X, <)=-|.e-*'. f ' (^)./, [^. \J 1?- (“)*] ‘dl. .... (308) 

X-^ 
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(«+«<)+/ {no-ut)] 


u-u . 

+ — -e >ci. 


r+wi 

•••• (309) 


Ul 


Into the expression for v {x, t) we shall now, according to (294), 
insert 

-f{i)^lc-p-F{i)+(al+rc)-F{§)+lc-G{§) .... (310) 

or, as follows from (299), 

-/(^)=i[p-^(^‘)+2«--P’(^)+G{^)] (310a) 

U> 

Hence formula (308) gets the appearance 

V {x, t)= 

x+tU 


'‘fi<l2‘ 




If we here compute the derivative according to (309), we arrive 
at the following ultimate expression for the voltage function: 

V (x, t) —^”^ 2 ' (x—ut)+F (*+«<)] 

[«-^(i?)+G(s‘)]-/{>l-\/<^-(^)“]-4^ (312) 


iT+m 
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In this form, but deduced in another way, the solution of the 
problem is presented in Webster’s textbook of partial differential 
equations. We observe that, besides the exponentially attenuated 
wave process given by the first term of (312), there are also 
present constituents of quite another nature. To the correspond- 
ing value of the function at any point x will contribute the 
whole set of values of the functions F' and G within the inter- 
val X ± uL Hence we see that the procedure will not be like a 
wave travelling over the point, the effect of which may set in 
at one moment and cease again after some time, but will con- 
tinue for ever after the moment of beginning. Also the latter 
terms of the expression are furnished with an exponential damp- 
ing factor with respect to lime. 


E. A disturbance problem of the telephone technic. 

As another application of a direct practical interest we will 
now treat a special case of the general transmission-line pro- 
blem outlined in Chapter 6. Thus we take under consideration 
a homogeneous telephone line of length equipped with some 
apparatus at the ends, the line being suddenly exposed to the 
sinusoidal magnetic field of a parallel power transmission line. 
The disturbance field is supposed to be the same along the tele- 
phone line; it will cause a voltage transient on this, the nature 
of which we are going to examine. 

The functions CP and 'F in (96) now get the very simple form 

<D {x; p)^0 {p)=-dp- [B {p)-BiO)]=J^^-dBo, 

p-jrci) 

V(x,-p)=0. 

Then we have assumed that B (t)=:Bo sm wt H (<). 

By (98) and (100) the voltage equation will be 

/i“F 

-±-dHp).F=0. 

In accordance with the methods of Chapter 6 we want a solu- 




no. 


tion of this homogeneous equation under the inhomogeneous 
boundary conditions 


where 


F(0,rt-E’.W+^.f(0;p): 


(314) 


e'Ap) 


E\ip)- 


Z' (p) 
r+pl 

Z''{p) (op^ 


r+pl p^+io^ 


'd‘ Bo 


•d*Bo. 


(314 a) 


No inhomogeneous differential equation occurs in this special 
case. 

The complete solution of the problem we have already got 
expressed by the formulas (104 a) or (125). For the sake of 
simplicity we will confine ourselves to the simple case of a dis- 
tortionless line, the ends being circuited by pure resistances. Then 

d=-{p+^0. Z{p)=R, Z'{p)=R', Z"(p)=R" .... (315) 
u 

and, according to (121), the coefficients of reflection are pure 
numbers, in magnitude < 1 ; 


„_R-R „ R-R" 

R+R' ’ ' R+R" 


(316) 


Further, the electromotive forces are Ip 


EVp)= 

E"*(P)- 


P 

P+P 
. V 


lop 

(op 


p-\-P 


■ I 
R" 
" I 


d* Bo 


* d * Bo* 


(3.17) 


The corresponding primary functions may be written down 
immediately on account of the formulas (74) and (75). Hence- 
we obtain, adding the factor H {t) to each side, 
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e'*{<)= •(— w<+wsin o)t)" d-Bo-H (t), 

(318) 

0 ) 

(0=" 'a ^~2 ■( — /‘-e^i^^+ZH-cos wi+w-sin -d-Bo-H (t). 


Now, by (125) 

v'{x,t)=W'{x;p)-e'^(t) and v" {x,t) = W" {^x; p)-e'\{t), .. (319) 


where W and W" are certain series expressions formed by terms 
of the type 


R 

R+R' 


G'fc-G''* . 0 — ^ • (2A:s f jc) 


Operating upon e'* (<)> this integrator gives: 


R 

R+R' 




2fc^+a? \ 
^ /' 


(320) 


Hence we obtain a system of waves, exponentially attenuated 
and travelling along the telephone line in both directions with 
the velocity u. A reduction of the amplitudes of the waves is 
also caused by reflection-losses at the ends. In practical appli- 
cations the waves are often damped so rapidly that even a few 
terms offer a sufficiently good picture of the transient disturbance 
process. The steady state disturbance is easily obtained by 
omitting the exponential term ii> the expressions for (f) and 
e"* («). 

The general case, where the line is not distortionless and 
where Z' and Z" are functions of p, can always be attacked 
formally in the same way, but we see at a glance that con- 
siderable complications will enter. 


F. Application of integrational and operational methods to the 
double-pole-radiation problem. 

The fundamental equations (Maxwell’s equations) of the electro- 
magnetic field may be written: 
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rot =47r l+^. 

at 


rot E 


dB 

' dt 


(321) 


ff, D==e E, /=a (E+B,,). 


(322) 


H denotes the vector of magnetic field strength, E that of 
electric field strength, B the magnetic induction, D the electric 
displacement, I the current density, and Ek the external electric 
field strength. f.i is the permeability, e the dielectric constant, and 
<j the conductivity. 

We shall assume Ek to be zero and e, and g to be constant 
all over the space. Further we prescribe the fields to be iden- 
tically zero at the moment /=0. Thus, if henceforth we let the 
capital letters represent the basic integrators, we get the integra- 
tional equations: 


rot H=(47ra+pe) 'E, 
rot 


(323) 


From this the solenoidal nature of the vector fields necessarily 
follows, i. e. 

div E=0, div (324) 


We will examine the process of radiation from an electric 
double-pole, the axis of which we take as one axis of a spheri- 
cal coordinate system (r, 9). From the cylindrical symmetry 

with respect to this axis follows that all magnitudes will be iii' 


dependent of the angle (p, and thus we may put ^ the 

equations. In the very neighbourhood of the doublet, or the 
))rational current elements) as it is sometimes called, the finite 
propagation- velocity of the disturbance can be let out of conside- 
ration, i. e. we may as a first approximation assume the fields 
to change simultaneously at all points and to coincide with the 
fields in the so-called quasi- stationary case. In this the compo- 
nents Ey), Hn and are all zero, which means that the elec- 
tric vector lines are situated in planes through the axis, and the 
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magnetic lines are circles with their centra upon the axis and 
their planes orthogonal to it. The magnitude of the vector 

H, H=H^, is connected with the current in the element 

according to Biot-SavarVs formula, which we may write in the 
integrational form 

H=~l p Q (325) 


Q being the basic integrator of the charge and I the distance 
between the two poles. 

We shall make tentatively the assumption that the general 
radiation field of the doublet has the Same fundamental geo- 
metrical character as the one just outlined, and especially that, 
when the distance r tends to zero, the analytic form of the mag- 
netic field approaches the expression (325). Then, according to 
well-known formulas as to the operations rot and div in spherical 
coordinates, we get from (323) and (324) the following system 
of partial differential equations: 


^ ~{sin^-H)=^{4na+pE)-Er, 


r sin^ 


-y^{r-H)={ina+pB)-E», 

1 IP \ 1 TT 

sin^~(r*'£r)+r- ^(sin ^■£#)=0. 

or o^ 


(326) 


From the last equation we conclude that Er and E» can be 
derived from an auxiliary function V. Indeed, it will be identi- 
cally satisHed if we put: 


Er-- 


1 


E»= 


du 


1 


dJJ^ 

r-sin -9 dr 


(327) 
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Then, from the second equation (326): 

. 1/- . m = ^ . 
T dir ' r-sin^ dr 


Here we may divide by - and then integrate both the members 

with respect to r, omitting the constant of integration (Cf. (326)). 
Thus 

(328) 

rsm^ ' 


Inserting the expressions found for H and Er into the first 
equation (326), one finds this to be identically satisfied. From 
the remaining equation of the system we get at last, putting 

(5^(p)=p*-fM+p-47r<7|tt (329) 


a partial differential equation for V, viz 


d^U sin ^ 

gj.2 + J.2 


dll du\ 

fil'Sysin & j 


-d“-C/=0. 


By the substitution 

17=sin=‘^-a)(r) 


(330) 

(331) 


one can eliminate the variable ^from this equation, thus obtaining 
an ordinary differential equation for <P 



• O. 


(332) 


The field components, according to (327), (328), (329), and (331), 
are now given by: 


0 

H=—-sin^~ 
Pitt r 

0 

jEr=2cos^*— , 


0 ' 

JS, 9 ^=— sin ^ — . 

r 


(333) 


Two independent particular integrals of (332) are known to be 



115 


0i=Av(^~+d\e-\ ©2=^2-/- —dVe^’’ 


The factor Ai being determined so that H satisfies the condition 
(325) at the origin, we find, corresponding to (Pi: 


H=sin 

Er=^-cosO-l~i+^^e-^’^-l-p-Q, (335) 

,9.(^+A+1J 


Replacing 6 by —6, we get the system corresponding to 0^- 
In the most important case the surrounding medium is an 
isolator, i. e. (7=0. Then, from (329), if we choose a suitable 
branch of the square root: 

<3=H „=-L. (330) 

u \Jen 

In this case we have the interpretation 

e-^^- p Q H {t)=e~^« q {t)=q'{t-^ (337' 

the right hand side representing a spherical wave travelling out 
from the doublet with the velocity u. The other branch of the 
root corresponds to a similar wave converging towards the ori- 
gin ; obviously such a wave cannot exist in our state of affairs. 
Accordingly, attaching the unit-function H (t) to each member 

and interpreting p as the differentiation-symbol—, and at last 

Ov 

denoting for short by m the moment of the doublet 

T 

at the instant we obtain the following explicit expressions 

for the field components : 
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TI • I* /W” , ^ 

£[=sin ./•( h 


jEr=COS ^ 


ru ’ r' 
'2iiiurn 


/ It 

t^'in 

f^u-m 



1 “ y.2 



.... (338) 


Because, according to our assumptions, g (^) ~ 0 when / < 0, 
there is no electric or magnetic field strength present at a point 

V 

at a distance r from the doublet, until the time — has elapsed. 


When the conductivity of the medium is not zero, the radia- 
tion problem becomes much more complicated. We shall now 
briefly outline how it is possible to attack this problem by inte- 
grational methods. One general way of treating a compound 
integrator is to split it into terms or factorials of simpler forms 
or of standard types the interpretation of which is already 
known. For that scope we first rewrite the expressions for the 
field components under the following form: 






Er=2l^-l-COS p .ff (t) 


jE^=/t-Z-sin^‘ 


r d i> 


(,)+LE.Q.E.^ . H (0 j . 


.. (339) 


Here all the factors appearing are familiar integrators. Because 
of the formulas (297), (299), (303), and (329) we have: 
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d{p)=--\f{p+Ky-A\ ^ (340) 

V/ 

1 , 27ra 

U=—p=r, K = /i = . 

ysfi « 

For the special value r=0 we get the interpretation of another 
V 

integrational factor, ^ . Hence 

(341) 

At last we meet in (339), besides Q (p), which we may regard 
as a given basic integrator, the factors and p, symbolizing 

pure differentiation processes. Indeed, from (329) follows that: 


It may be sufficient to carry through the method upon one 
of the terms. We choose one of the most complicated, say the 
second term of the expression for the magnetic field strength: 

(m 

If Q and*F denote the integral operators (?(0)=0) 

t t 

{t-r)-dT], ‘f^=[y'(0)+yV'(<-r)-dr], .... (344) 


this can be presented in the more explicit form 


i-sini9 / d 


• ^ + 47ra/t ) • X Q X Z (O' 



(346) 
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From (341) we see that y){0)—u, and for t>0 

Thus we know the expression for the operator W. 
If we put v(t) = Wx x(t), we get: 


H* 


I sin -9 


£jm-9'(0)u(<)+ 


t t 

+ 47r(T|if J"q\t~t)'V {r) dr+€|it' f{t—T)'V{r) drj, 
0 0 

t 


(347) 


-{~X'T 




Accordingly we get a wave with the propagation-velocity u and 

r 

— X 

a spatial damping-factor ^ ^ of the face. The function v (t) 

characterizes especially the mechanism of radiation. 


12. Some applications of integrational methods to linear 
integral equations. 

In Chapter 6, when studying the universal equation of tele- 
graphy, we arrived at a differential equation of the type: 

dW 

j^-d%p)V=—F{x;p) (348) 


In order to find the solution of this equation corresponding 
to the homogeneous boundary conditions 


'"(«> P) r+pl dx r+pl dx 


{s;p).... (349) 


we first tried to determine the so-called Green function K{x^§;p) 
belonging to the homogeneous equation under the same bound- 



119 


ary conditions; having found this, we could form the solution 
desired according to the formula 


V (a>; P)=J'k {CO, 1 / v) F {§; p) d§. 


(350) 


Another mode of proceeding will bring us into the domain of 
linear integral equations. By this we take under consideration 
the simpler differential equation 

^^-U{x;p) (351) 


and try to solve it under the same boundary conditions as be- 
fore. Denoting by H {x, <?; p) the Green function of the differential 


expression 


we can write the solution of (351) 


V{x)^ /Hix,i;p).U{§;p) dl 


(352) 


We can apply this result to the original equation (348) bj^ 
introducing into (351): 

U {x; p)=F {x; p)—d%p) V{x). 

Then the equality (352)* will appear in the form of a linear 
integral equation of Fredholm's type: 


where 


s 

V{x)+d^{p)- J'H{x,§;p)-V{§)d§^R{x;p) (353) 

0 

8 

B{x;p)=J'H{x,i;p) F{i;p) di. (354) 


The kernel of this equation, H {x, §; p), is symmetric in x and 
but it is real only provided that p is real. One finds without 
difficulty the expressions: 





«•/„ JZ' {p)+{r+pl)’a!] [Z" {p)+{r+pl)-(s-§)] 

a I®. V) ^2! (p)+Z" {p)+(r+pl)-sHr+pl) ’ 

if !e<i 

[Z-(p)+(r+pi)-s1-[Z>)+(r+ 

^ .c, V) {v)+(T+viysY{r+vl) ’ 

if 


. . ( 355 ) 


Therefore, if Z' (p) and Z" (p) are rational functions of p, this 
kernel is itself such a function, and accordingly it is of a simpler 
structure than the Green function K {x, §; p), which involves 
transcendentals in p (Cf. (107)). 

The treatment of the equation (353) may be executed in 
accordance with the general methods of the theory of integral 
equations. From the point of integrational calculus special com- 
plications will enter, because the kernel H {x^ §; p) contains the 
parameter p. In the following we will discuss more in detail a 
special case of practical importance where that is not so. We 
arrive at this case putting 

Z' (p)=5'.(r+p0, Z" (p)=5".(r+p0, 

s and s' being positive numbers of the dimension of a length. 
A transmission-line can always be so connected at its ends that 
this condition is fulfilled. Then we have simply 




^ -i-x) {S -ts- 
/+s'+s 

if x^^ 


H{x,§)- 


s'+s'+s 
if x'>§ 


Now the integral equation (353) has a real symmetric kernel, 
and we know that the corresponding homogeneous equation 


(p {x)=X-J'u {x, §)-(p {§) d§ 


(357) 


has real characteristic numbers^ which may be proved to be all 
simple and positive. We will denote them by and 

the corresponding normal junctions by 99i(«r), {x) 
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According to the so-called Schmidt's formula we get the following 
expression for V{x): 

oo 

V (358) 

n=i 


putting : 


■Rn ® 


p)Mi)di (359) 


By the definition formula (354) we obtain, making use of the 
symmetry of H {x, i) and of (357), 


S 8 

Rn ip) V) 




(360) 


Thus Rn{p) differs only by the factor from the Fourier- 

71 

coefficient of the function FJ^x; p) with respect to the system of 
normal functions (pn{x). 

In our case we can easily • get some knowledge of the 
characteristic numbers and the normal functions. Indeed, put- 
ting F {x,p)=0y we infer that the integral equation (357) is 


equivalent to the differential equation (^=(o®) 

(3ei> 

combined with the boundary conditions in question, viz 

(p (0)=5'*9)' (0), —(p{s)=^s''(p{s) (362) 

Writing the solution of (361) 

(p (ic) =A • sin CO x+B • cos w x. 
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we get, according to (362), the following system to determine 
the coefficients A and B: 

(sin ws+ws" cos ws)-^+(cos cos — ws'' sin cos) B=0 
cos'-A+{—l)-B=0. 


The determinantal equation of this system, i. e. 


tg tos 


w(s'+s") 

1 ’ 


(363) 


gives us the possible values of (o or A. 

The equation being transcendental, the properties of its roots 
may be most easily overlooked by a geometrical discussion. It 
means no serious limitation, principally at least, if for that pur- 
pose we assume s=s'—s. Putting ws=x, our equation becomes: 


tgaj= 


2a? 


(363 a) 


y 



Fig. 10. 
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Referring to fig. 10, we note that A=ci)^=— takes a sequence 

s 

of positive values which asymptotically tend to the form 


An = OJ^ CO 


7t^ 


(364) 


when n tends to infinity. Hence the normal functions are 

(pn{x)=An‘{sin OJnX+Ci)n-S -cos COnX), (365) 

and An may be chosen so that cpn {x) is normalized. To each 
characteristic number one normal function corresponds. The 
orthogonality of the normal functions can be verified e. g. from 
the symmetry of the function H («r, <?). 

Now we return to the expression (358) for V (x). Attaching 
the unit-function to both the members we obtain the time- and 
space-function v {x; t) in the integrational form 


v{x,t)=R{x;p)-H {t), . . (366) 

n=i 

which shall now be interpreted. 

At first we discuss the case where no electric or magnetic 
field acts upon the line, but there is originally a certain distri- 
bution of voltage V {x, 0). Then, according to the definition 
formulas (100) an (96), 


F {x; p)={p^'lc+p rc) v {x,0) (367) 

and, according to (360), 

8 

Rn (p) = ip'^ lc+p rc )-^^ ; 0) (pn {i) di 

» .... (368) 

-ip^ lc+p rc) 

U 


Further we obtain from (354): 


s 

R (x; p) H {t)={p^ lc+p rc)- ^ S {x, i)-v {i, 0) d§ H (0=0 


(369) 


as, for t>0, p is equivalent to 


dt 
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Therefore 




In this expression appears a function 


(371) 


which obviously may be regarded as characteristic of the pro- 
blem. In the universal case we have (Cf (299)) 


1 




u- 


and thus 


if we put 


70m (<) = |l— 


('P + uf + iO, 


■H{t) .... (371a) 




(372) 


Employing the expansion theorem, we readily arrive at the iden- 
tity: 

(p+«)2+-,.--ff(0= 

(373) 

1 r __ _ -1 

| ^ 7^ — (wn-COS (Ont + K Sin WnO J -ff (t), 


co: 


from which for positive t we get the derivative 


col 


(I)n 

Accordingly, as !x?n(0) = l, the integral operator corresponding to 
Wnit) must be: 

t 

jr„=|^l-7i"-^yV-*«-^)-sin w„(<-r) drj (374) 


We will soon make use of this result. 
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The general term of (370) involves the integrator 




{p+Kf—A’‘ 

(P+k)"+< 



(Cf. the definitions of A and n in formula (299)), the primary 
function of which may be obtained from (373) or directly by 
the expansion theorem. Thus we get from (366) the result 


oo 

V (x, t)=^Vn-(pn (x)-e~’‘*-^OS <Ont+^- sin (Ont^-H (t) . . (375) 

w=l 


If we define an angle cc according to the equations 
A==U'OJn‘COS a, (On=u ajn' sin a, 
our expression for the voltage can also be written: 


V (x, 0 = ^' 


= yu—-Vn-9>n(x)e~‘’‘-sin(wni+a)-H(t) (376) 

On 


When t—O, as is most easily seen from (375), this expression 
is reduced to; 


V{X, 0)=^Vn-(fn(x), 


(377) 


which is the Fourier-development of v (x^ 0) with respect to the 
system of normal functions (pi{x), (p 2 (x)y 


Now we will pass over to a brief examination of the case 
where the voltage as well as the current is initially zero, while the 
line at the moment is exposed to e. g. a magnetic distur- 
bance. Then, according to (100) and (96), 

F {x; p)=d p [B'x (x; p)-^B'x {x, 0)] (378) 

We will limit ourselves to the practically most common case, 
where the external magnetic field can be regarded as varying 
simultaneously along the line. Then we may put 

B {x, t)=B {x)-a (^)==B {x)-A {p)'H (t) . 


(379) 
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and 

F(a;;p)=d-p-B'(a;) [^(p)-a(O)] (380) 

Further we have by (360) 

^n{p)=^ d-p-[A (p)— a (0)] B'n (381) 

n 

employing the notation 


B'„= / (382) 


and from (354) follows: 


S 

R (cc; p)=d p [A (p)-a (0)]- J ' ff {x, <f) B' {i)-d§. .... (383) 


Thus we obtain, making both members operate on H (t), 


Rn{p)H{t)=-^„-'B'na{t), 


s 

R (x; p) H {t)=d J'H {x, ^) B' {g)-d§-a (t) (385) 

0 

Inserting these results into (366), we get the following expres- 
sion for the voltage: 

V {x, t) = 

g CO (386) 

d J'H {X, g)-B' {i)-dS-a' it) - it)> 


and can now immediately apply our formula (374) to arrive at 
the wanted explicit form of the function v{x,t)> which becomes: 
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n=i 


V {x, t) = 

t (387) 

(pn (x) • ■ sin Mn (t—r) a' (r) dr. 


as we have, because of the integral form (354) of B (x; p), 
8 oo 

J~ {X, i) • B' (^‘).d.^ <pn {X). 

0 n=i 


In passing we remark that the function v {x, t) derived above 
only represents one constituent of the real voltage transient. In 
fact, because of formula (103) we have another constituent ob- 
tained from the homogeneous differential equation under inhomo- 
geneous boundary conditions. The examination of this part does 
not lead to integral equations and may be performed according 
to methods exemplified in the preceding chapter. 


As a last example of the treatment by integrational methods 
of a problem where initial and boundary conditions are combi- 
ned, we shall derive the analytic expression for the motion of 
an inhomogeneous elastic string with fixed ends. 

The differential equation for the elongation v (x, t) will then 
be of the form • 



where a [x) denotes the modulus of elasticity multiplied by the 
section, a {x) the mass pro unit length, and e (Xy t) some func- 
tion of disturbance. 

Our mode of reasoning may as well be applied to the more 
general equation 

{x) v=a (*)-^+e {x,t) (389) 

Here the left hand side represents the most general selj-adjung- 
ated linear differential form of the second order. The functions 
a {x)y a {x)y and (x) are prescribed to be continuous in the in- 
terval 0<x<s. We assume the initial conditions 
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v{x,0)=F {x)y ~ (a?, 0)=G (iU) (390) 

and the boundary conditions 

i;(0, /)=0, v{s,t)—0 (391) 


Then, of course, we must have jP (0)=jF (5) = 0. 

From a purely methodical point of view we could as well 
have chosen other homogeneous boundary conditions. 

Now we introduce the basic integrators of v (x, t) and e {x^ t)^ 
putting 

V {x,t)=^V {x; p)'H (t), e {x,t)=E {x; p)-H {t) (392) 

and accordingly obtain for V (x; p) the ordinary differential 
equation 

^(« (x)V-p--a i^yv^r {x; V) (393) 

where 

y (x; p)=E (x; p)—p^-a (x)-F (x)—p-a {x) G (x). .. (394) 

Further, let G {x, §) be the Green function of the differential 
form 

(395) 

under the boundary conditions in question. 

Performing the same reasoning as above when studying the 
equation of telegraphy, we easily prove the following integral 
equation to be satisfied by the function V {x;p): 

s 

V{x)+p‘-J'G{x,i)a{§) V(i)d^-=q{x;py .... (396) 
0 

where we have put for brevity 

s 

(l{x;p)=—J'G(x,g)-y{§;p)-d§. (397) 

0 

Unlike (353) this equation has a kernel G{x,i)-a{§) that is 
not symmetric, unless the density a {§) is a constant. However, 
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one can readily pass over to an integral equation with a symmetric 
kernel also in t his case. For that scope we only have to multiply 
(396) by {x) {a {x) is supposed not to assume negative values 
in the interval (0, s)) and then to introduce the new functions: 


V^(x) = \fa{x)-V {x), 

G* ^ = |/a (aj) Ya (|) • G (x, 

Q* P)=ya {x) Q {x; p). 

Then our equation will appear in the form 


(398) 


V^{x)+p^^- / G^{x,Syj\{iyd^=Q^{x;p) (399) 


The kernel now being real and symmetric, because the Green 
function G {x^ has these properties, we know that the corre- 
sponding homogeneous integral equation 


(p{x) 




G* {x,§)<p{i)d§ 


(400) 


has an infinite set of real characteristic numbers. These may 
be proved to be all simple and positive; we will denote them 

by o)\, col corresponding normal functions by 

<fh{x), (pi{x) 

According to Schmidt’s formula the solution of the inhomo- 
geneous equation (396) can be written: 

oo 

V* (a!)=Q* {x; {p)<pn (x) (401) 

W-1 

where 


8 



0 


(402) 
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Now we make both the members of (401) operate upon the 
unit-functii^i H (<). From (394), (397), and (398) follows: 


8 

{x; p) = j/'a ^ («. ^) [?>®-a (^)--P’ (^)+ 


, (403) 


+p-a(^)G(^‘)-i:(if;p)]-d^‘, 

d 

and we get, because P—^ t>0, 

8 

Q* (x; p) H (<) = -)/^) J^Gix, i)-e {§, t) di (404) 


Further, by (402) we have 

8 S 

Q»(P)= /Iv/^)- f ^ii,r;Hp^-air,)-F(r])+ 

» “ (405) 

+p a (»/) G(»?)— E (ij; p)] d»?j-9?n 
Hence we may write: 

Qn(p) = p'''Fn+p-GK—En(p) (406) 

Employing the notations 

F:,(x) = y^)-F(x), G*(x)^^^) G(x), (407) 

we can define the coefficients Fn and Gn by the double integrals: 
8 8 

^»=yy' SJ* (<f, J2).F* (rj).^n {i)-d§ dri, 

0 0 

8 8 

f V)-G* inYVr. {§)-di-d>j. 


( 408 ) 
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On account of the symmetry relation and 

because (pn{§) satisfies the homogeneous equation, we can per- 
form one integration, thus obtaining the simpler expressions: 



The last term of (406) is defined by the integral 

8 8 

En{v)=J'J' G {n; p) -^n{^) dldr,. .. (410) 

0 0 

With our new notations the formula (401) presents us the 
function v (x, t) in the form 


V (x, /) = 


8 

—sjalxyj'i 


G(x,i)-e(S, t) d§ 


-^(p^.Fn+pGn)-:^,-H{t)-<pn{x) 

n»=i 


. .. (411) 


n==i 


At last, making use of rule nr 5 of our table and the nota- 
tions (409), we arrive at the explicit function v (x, t): 
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8 

V {x, t)= -\la{x)-J' 1 {x, ^) e (?, i) d^ 


+ N (i»)*COS (Ont'H + -(pn(0[!)- sin (Ont-H {t) (412) 


n==i 

t 


+^9>n(®) 

n=i 


) — iOnJ ' ii 


en{t)—iOn- / sill CDn{t—r) en{T) dT 


Here the first term may be developed into a series of normal 
functions, thus being found to cancel the first part of the third 
sum in (412). Therefore, on introducing 


{t) ■ 


A (»;) 




(t) 


(413) 


which expression is easily obtained by applying En (p) to H {t), 
we get the final result: 


oo I 

V {x, t) = — (®) 'J ' (^) 


(414) 


+ 7J^’n'«>n(a?)'COS (i)nt-H{t)-i-y ”' -(pn{x)-sin Ont-H (t) 


n=i 


n=i 


This involves, as we see, the Fourier-coclTicients of the functions 

}^a (x)-F{x), |/a {x) G {x), and ^ 

Ya{x) 

which, together with the boundary conditions (391), specify the 
problem. 



Appendix. 

1. The Fourier-Mellin integral theorem. 


Assume the real function / (t) of the real variable t to have 
the following properties: 

1. f (t) and its first derivative are piecewise continuous in 
each finite interval. 


2. / (0 “of/ — ^)] ^ point of discontinuity. 

mi 


3. The integral 


-TOO 


dr exists. 


Then the classical Fourier theorem states that the function 
/ (/) can be identically represented by the double integral 


30 + OO 

f (r) cos 

0 — CO 


( 1 ) 


which may also be written 

-}- oo 


-f oo 




— OO OO 


(la) 


provided that the right hand side of (1) converges if the cos- 
function is replaced by the sin-function; its value then equals 
zero, because sin is an odd function. 

By introducing the new variable p = i7, we can also express the 
contents of (la) by the two formulas 

-f lOO 

— too 


( 2 ) 
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+ 00 


F(p)=p T) e-^^ dr, 


(2 a) 


the path of integration in the first integral being the imaginary 
axis of the complex p-plane. 

As to the proof of Fourier’s integral theorem we refer to text- 
books of the subject, e. g. Courant- Hilbert, Methoden der ma- 
tematischen Physik. 

The Mellin modification of this theorem, as formulated in 
Chapter 1, formulas (1) and (2), is now readily obtained in the 
following way. 

In the right hand member of the first of these formulas, which 
by the second may be written 


a-\-‘ioo -\-oo 

R. M, (0 dr, 


we put p~o~\-iA, a being a real constant, A a real variable. 
Hence dp==i dA, and our path of integration runs parallel to the 
imaginary axis. Thus we obtain 


-j- CXD -f- OO 




(t), according to Fourier’s theorem (la), provi- 
ded that the function / (^) e~*^^ ' fulfils the condition 3 of inte- 


+ 00 


grability, i. e. that^^|/(T) exists. When this is the 


case 


— OO 


if ai < a < tta, any line parallel to the imaginary axis in the 
band region ai <Re{p}<a .2 may be used as path of integration. 


A theorem, reciprocal to that just proved, is also due to Mellin. 
This may be stated as follows: 

Let F (p) be a function of the complex variable p such that : 
{p = o+il) 


Fip) , 


P 


is holomorphical in the band region a<a</J. 


1 . 
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dl converges in the same region. 


3. I p I uniformly, when A ± oo in the region 
a+(i<a </? — dy (d>0). 

Then formula 2 of Chapter 1 is a consequence of formula 1. 
We refer to the book mentioned above for a more complete 
examination of these so-called reciprocal theorems. 


As a rule our investigations are to be performed on functions 
/ {t) that have no other point of discontinuity than possibly the 
origin; for ^<0 they will vanish identically and for < > 0 be 
more or less regular. Naturally, the integral theorem ought to 
be proved comparatively simply for such functions. We will 
now derive the theorem in a form suitable in many cases, and 
in a way that will give us a wider view of the notions of basic 
integrator and integral operator of Chapters 2 and 3. 

We start from the expression % 



(p{t) being any integrable function and j [t) a function having 
continuous derivatives in the Interval of variability employed. 
We are going to perform a suite of integrations by parts, and 
for that aim we introduce the auxiliary notations 
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As all the «n(0) vanish, the repeated integration by parts gives 

«/'(<) = 

= / (0)-Cp{t)+f (O)-Mi (<)+.... +/<"-*> (0)-«„_a (0 + 

t 

+ J'f”\t—T)Un—x (f) • dr, 

0 

and therefore we gel, according to (5), 

W (<)=/(0)-y'W 

+ J' + + 

0 

(6) 

t T 


The last term may be transformed in the above manner, thus 
getting ^the form: 

i t 

\^f \f j- 9 ’(r)dT (7) 

~~ 0 r 

Obviously, the function xp {t) can be regarded as derived from 
cp{t) by a symbolic multiplication (Cf. Chapter 3): 

ip{t)=:F'X(p{t). ( 8 ) 

the integral operator F being 

t 

^’=[/(o)+ y'{/'(o)+r(o).'-f^^+.... + 


t 

T 


(9) 
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The case n—\ corresponds lo the original formula (3), which 
gives 

t 

0 

i. e. the formula (36) of Chapter 3. 

Now we return to formula (6), specializing cp (t) lo the unit- 
function II {t). Then we know d priori that yf (/) i— / (/), provided 
that f (t) = 0 when / < 0. All integrations but one may now be 
performed and we get the simpler formula 

/( i) = {/ (0) +/' (0) ■ + • • • • +/<'-'>(0) - it) 

( 10 ) 

t 


Now, it* is possible to prove by means of the theory of residues 
and Jordan’s lemma that, if m is a positive integer or zero and 
a>0, 

(j+ioo 

fm 1 Z' 1 pvt 

J \;/-~’dp ( 11 ) 

\m 27r/ J p"'* V 


Accordingly we obtain 


u\ 


I (0)+'’,™+ ■ + ■ ■ (12) 


The rest-term may l)e transformed in the following way: 

t 

I (<- O-dr 


(t-r)dT- 


0 
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a+ioo 




(T + 'ioO oo 


-i/I/' 

a — ioo 0 


/^(ir).^ - dp. 

r 


pPt 


(13) 


The last operation, viz the inversion of the order of integration, 
is correct, if the double integral converges absolutely.* This will 
certainly be the case, if belongs to class nr n, and if a>a 

(Cf. the definition in Chapter 1; when In- 

deed, putting p — a + i • co, we get as an upper bound of the 
magnitude of the last expression: 



drl- 


d(i) 




Even the outer integral must converge, if n>l. 

If (13) be introduced into (12), we get the result: 


<f+ioo 



u 


This equality constitutes the integral theorem in the wanted 
form. By our deductions we have got it proved for any func- 
tion of class 2, but as we know, it also subsists for functions 
of class 1. The last equality (14) is, of course, nothing but 
formula (20) in Chapter 1. 


By a quite analogous proceeding we may generally prove the 
identity (Cf. Chapter 3) 


FXcp{t)^F{p).q>{t) 

* This condition is sufficient but not necessary. 


( 15 ) 
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to be valid for any function cp (t), the basic integrator of which 
(Cf. Chapter 2) may be In fact, we only need to employ 

the rule (Cf. page 36) 


n 




(t — (p{T) dr=~- 


(T+icx> 

/®(p) 


27ii 

a — ioo 


pU p 


dp, (16) 


which is a direct conse([uence of the rule of integration on 
page 15. 

The expression ip{t)= FXcp(t) has the meaning defined by (6). 
From (16) follows that 




cr-i-ioo 


v\ 


m , /"(o) 


/(«-)) (0) 


2711 

<J — J CO 




V 


pVt 


(17) 


where, in accordance with (16), we may write 



(T+ioo 

* 

a —ioo 


Transforming as above, we get 




K.= y dp^d. 


+ or; 




pPt 


{r)e~“'dT\<D(p)-^^dp. 


a — I’oo 0 


(19) 


Because | (D (p) | < A when o) > U, the absolute convergence 
of the double integral can be proved just as above, if ^i>l. 
Adding this result to (17), we get, according to (14), 

a+ioo 

V^(0=^y' ^iP)-^{p)ydp. 

(T -ioo 


(20) 
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But here, by deliiiilion, the right member agrees with the 
function F {p)'(p {t)=F {p)'0 {p)'H (t). 

The expressions (9) of F and (14) of F (p) may conveniently 
be termed the transformed operator resp. integrator of order n. 


2. A proof of Jordan’s lemma. 


Let O (p) be a function that fulfils Ihe conditions slated on 
page 10 in the region Re (p) > a. The path of integration may 
be a half-circle C of radius q and with its centre at m, where 
Re (m) > a. Then upon C we can put 


Denoting by e an arbitrarily small positive number, we can 
always find, on account of our assumptions, a number qo (f) 
such that the inequality 


I ©(/,)! <^ 


holds at all points upon C, as far as () > qo. Then we get the 
following estimation: 

3>(p) e^^ dp I 
c 


A I 


0 \ i Hin 6) • t 


Q’dO 


or: 


I/' 


0(p)’e^^'dp 


< 


2f 


0 

* ^ 7?e ( -t / « • cos 6 • t 7 /> 


TV 

Now we introduce 0=9'—^^ and make use of the fact that 

2 7t 

sin0'> — *0' in the interval O<0'<-. Hence we obtain, if /< 0: 

TT 2 


•2 

(l>{p)e^Hp ' ^ dff- 


2 Ue[m}.t jj. r V.jrl'a 

'^V" Jo’ 


7t 


•£• qe 
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or 


I/' 


® {p)e^dp 


<-i.e 




( 1 ) 


Also this inequality subsists for all ^ (e). However small 

the value of \t\=^ 0 is, we can always choose an 8 so that the 
right hand side is < fi being an arbitrarily small positive 
number. In fact, if 

8 < 8i‘\ t\ e ' ^ , (2) 

the absolute value of the integral must be less than fi, and this 
for all half-circles with () > qo (fi, t). Obviously this result involves 
the statement of Jordan’s lemma. From (1) follows that the 
convergence is uniform in the interval t< 0, provided that 
Re[m) > 0. 


3. A proof of the theorems of page 8. 

These theorems refer to a function of the complex variable p, 
presented in the form of a definite integral, viz 


F{p)_ 


f {r) e-P^ dr, 


( 1 ) 


/(r) being a real function of class nr 1 at least (Cf. p. 7). 

• F (p) 

The first theorem states that — is holomorphical in the 
region Re{p)>a, 

To prove this we will first verify that the infinite integrals 

oo oo 

J f dr and J " r / (2) 

0 0 

are uniformly convergent if Re{p] > a© > a. 

This may be done by estimating the rest-integrals h' and /o' 
with the limits t and oo. We obtain: 


= ‘"'■dr 

t 


oo 

< / «)• {(To— /3\ 

J Go— a 
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provided that t >t* (Cf. (3) in Chapter 1). Here the right hand 
side can be made less than any positive number £ if ^ is made 
sufficiently large, say if / > >^ 4 !, U being the same for any p 

of the region Re{p}>ao >a. 

Further we have, if t >t*. 






<<ro— a)*T. 


{Oo-af ^ 

Obviously, as X'e~'^ tends to zero if x tends to infinity, the 
right hand side of (4) will remain less than e \{ t >12 U 
being the same number all over the region Re^p) > ao >a. 


Now, as h converges uniformly and its integrand r*/(r) e^^ 
is a continuous function of t and p, we may conclude that I 2 , 
but for the sign, represents the derivative of h with respect to p. 
F(p) 

Hence the function — — has a derivative at each point, and 
P 

accordingly it must be analytic anywhere in the region Re{p) >a. 


F(p) 

The second theorem states that tends uniformly to zero 

if I p I -► c>o in the half-plane Re{p} > Oo >a. For the proof we 
will employ the formula (14) in Chapter 2: 


F(p)==f{0)+J'f(r)-^^dT. 


Of course, we have only to verify that the last term of this 
expression is uniformly limited within the region in question. 
We obtain: 


00 



(T)e~^^dT 


0 r* 
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I. oo 

< M- j'e-^'? dx 


denoting by M an upper bound of | /' ( t ) | in the interval (0— t *). 
Accordingly 


\f< 




< 


M 1 

— •(!— (5) 

ao ao— a 


and the right hand side represents an upper bound of the mag- 
nitude of the integral, which is the same for all points of the 
region Re{p} > Oo > a. 


4. Note about the convergence of the modified integral. 

It is almost trivial to prove the uniform convergence for^ > U 
{to an arbitrarily small positive number) of the complex integral 



LiL^i 


0 (p) being a function that increases at most as a power of p, 
when p->oo along Li and L 2 (Cf. page 17 and fig. 2). 

Let us split the path of integration into three parts L'l, 

L'^i by the points Pi and P 2 , the first and the last part having 
an infinite extension. To get the uniform convergence established 
we only have to verify that it is always possible to fix the 
positions of Pi and P. so that the contributions to the integral 
from L'l and L '2 remain less in magnitude than an arbitrarily 
small positive number for all values of ^ > <0 > 0. Obviously, it 
is sufficient to perform the proof for one of the parts, say L'o. 

Upon L 2 we may put (a > a). Accordingly dp = 

dQ • and 

^pl _ • COS yj+iji • sin y) • t 

Therefore we get 

00 

1A1<^- J'\0{p)\ e9--o»'r-t.do, (2) 

(>9 

denoting by the contribution to I from L' 2 . 
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Because of our condition that | ® (p) | should increase at most 
as rapidly as a power of [ p | when ^ oo, there must surely exist 
a number such that the inequality 

I® (p) !<«*■•? (3) 

subsists for all values of ^ > ^» 2 *, k being any fixed positive 
number. Hence we obtain, taking p.. > 

oo 




r 

f'‘> l< .. / • P do (4) 

^nj ' 


This integral converges if k+t coscp <0, i. e. if 

t {— cos >k, 

and we get 

^ ^A;()u+((r-f (12 •cosf)»t j • cos (/) •to 

t{~cos(p) — k 2iz to{~ cos q)) ~ k' 


... (5) 


provided that 


t > to > — ^ — and 02 > - — . 

— cos (p — cos cp 


(R) 


Accordingly, we have got | /' 2 1 < an expression that may be 
made arbitrarily small if we take ^2 sufficiently great, and this 
for all values of t > to. It should he observed that the bound 
to can be regarded as arbitrarily small, since k could be any 
positive number. 

The integral I'l can be treated in exactly the same way, and 
hence the uniform convergence of I is proved. 

Fiv) 

In the modified integral (j?) on page 17 we have ®(p)=Q(p)- — 


where Q (p) is assumed to increase at most as a power of p and 
F(v) 

— — to tend to zero when p-*-oo in the region in question. 


Thus the uniform convergence of the integral (/?) follows a fortiori 
from that of I. 
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In order that F (p) shall be uniquely defined we assume cuts 
e. g. in the form of straight lines to be put between the branch- 
points qi, qi -- and the point at infinity in the domain to the 
left of L\ and La; obviously these cuts can always be drawn so 
that they do not intersect. 

Now we shall prove the expansion theorem: 

oo oo 

/ tf)~F {p)H (‘) ( 2 ) 

1 V= 1 

At first we may note that, according to our assumptions, /(t) = 0 
when t<0. To evaluate the integral 


'<‘>-4/*’ '*'>■ 7 ''’’ 

LiJU 

when / > 0 we introduce: 

n 

^’(P) = 2^+/?»(p), (4) 


where, e being an arbitrarily small positive number, we shall have 

I (p) I < f 

for all values of n sufficiently large, say n > N (e), N being the 
same for all points of Li and L 2 . " 

Hence we get (^ > 0) : 


n 



L1L3 LiLi 


A rough estimation of the magnitude of the last term in (5), 
sufficiently good for our purpose, can easily be obtained, if we 
observe that on the path of integration we may put 

IpI > |(x| sin 9 ), \d'p\^dq 

Thus: 

LiU 




TT sin qj cosq) < |a| 


( 6 ) 
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Accordingly, t being any fixed positive number, the right hand 
side of (6) can be made as small as we like by taking e suffi- 
ciently small, i. e. by taking n sufficiently large. Therefore we 
have proved for / > 0 the equality: 


oo 



V— 1 LiTjj 


*Here the origin is the common branch-point or pole of all 
the terms, and obviously we may replace the path L 1 L 2 by a 
loop enclosing the negative real axis. Thus (2) follows imme- 
diately. 

We see at once that our previous rule 

/(0) = F(oo) 

for functions F {p) that are regular at infinity still subsists for 
the general class of basic integrators under consideration. 


Now we shall proceed to investigate more generally the 
asymptotic form of a function f {t)=F {p)‘H (t). For that aim 
we will first examine the contribution of a singularity at the 
origin by the same method as in Chapter 9, but this time assuming 
the cut to be any straight line from the origin to infinity in the 



left half-plane. (Cf. Fig. 12). Accordingly we take under consi- 
deration the integral 





(7) 
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extended over a loop S enclosing the cut. About the origin, say 
within the circle \p\<a, we prescribe F (p) to be of the form 


F(P)=^ 


ayjfy+p^-Rn 


( 8 ) 


where Si < Sr-- < s», s > 0, and such that 


|ie™|<e 


TtS 

a® 


(9) 


for any sufficiently great value of n. Further we assume that it 
shall be possible to determine two real numbers If > 0 and k 
so that the inequality 

(10) 


subsists anywhere upon the borders of the cut outside the points. 

Our mode of reasoning will now be identically the same as 
in Chapter 9 ; we only have to change t into t • cos cp in the 
formulas (191) and (199). Hence, by chosing n sufficiently great 
we must get 



•< Sf 


(11) 


for all values of t>ti and U- 

From the definition formula of a basic integrator we further 
conclude that 


(TH'ioo 



provided that Sv<l. If ^v>l we may regard the integral along 
S as the modified integral, making it possible to maintain the 
last equality for all real values of Sy, and thus it follows that 
(11) is equivalent to the asymptotic expansion (201) or (202) 
for I {t). 
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jP (v) 

The general case, where - - has an arbitrary but finite num- 
ber of singular points, can now be treated without difficulty. 
Obviously it means no essential limitation if, for the sake of 

F(v) 

simplicity, we confine ourselves to e. g. the case where has 

'P 

the poles po=0, pi, and the branch-points gi, q^. We also 
\F(v)\ 

assume I | to approach zero uniformly when p-^oo, so that 

Jordan’s lemma may be applied. Then we know that /(t)EEO 
when t < 0, and 


and 


/ (0 = 7,,,+/5.^+/p(,+/j>,+73J2 

(13) 

putting e. g. 

1 r 



Hi 

(14) 

1 r 



o Pi 

(15) 


Here Si denotes a loop surrounding a straight-line cut from 
the point qi to the point at infinity in the left half-plane (Cf. 
fig. 11). 

To compute /g, we change the variable of integration, putting 
p—z+qi, thus getting 








+9')- 


z j 

dz. 

5+9i z 


(16) 




The path of integration in the ;:-plane becomes a loop sur- 
rounding the origin, obtained by a parallel displacement of Si 
a stretch given by the complex vector — ji. Upon this path we 
may employ the results obtained above. For that aim we assume 
the e([uality : 

n 

F(z+q,)-^^^-=^a,.-z^y+z^-Rn (17) 


to subsist within a circle of radius a about the origin, the 
factor Bn satisfying the condition (9). Of course, we may as well 
prescribe that 
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F{p) 


■^■(p—qi)=^ar (p—qiY’'+(p—qiY Rn (17a) 

)-=. 1 

in the region |p— < a about qi in the p-plane. This being 
the case, we get from (16), as we know, the asymptotic expansion 






(18) 


or, symbolically written. 




^jir-p^y -H (t). 


(18a) 


One can arrive at this expression in a simpler way by employ- 
ing formally a shifting rule of our table in Chapter 5, viz : 


F ip) H {t)=e<‘*-F(p+a) e-<^i H (t). 


Here we put a=qi and 
getting 


. P_ 

p+qi 


H{t) 


instead 


of 


F{p).H {t)=e'‘^^-F(p+qr)-^^^-H (t) 


thus 


(19) 


After that we only have to employ the equality (17), replacing 
z by p, to recover (18 a). 

In a quite similar manner the constituent Iq.^^ corresponding to 
the other branch-point 52 , is obtained. 

The contributions Ip^, Ip^ from the poles can, of course, 
be computed by the same method, but in this case exact instead 
of asymptotic expressions are obtained. It may be of some 
interest to compare this mode of proceeding with other methods 
of computing Ip^ = (F {p)'H (t))p^. 

We assume F(p) to be of the form 


(P(p))pi==br-^-^^+bi-~^+ch-p+aip{p—pi)+ 
+afp (p— pi)“H 

in the neighbourhood of the pole. Accordingly (Cf. (17 a)), 
— -(p— pi)=52--- _|_p^ +5i+ao(p— pi)+ai(p— pi)H — , 


( 20 ) 
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and we get by formula (18 a) 

A . . . j .H(f) =ePx<. {bi+b, t)-H(t). (21 ) 

Alternatively we may apply the shifting rule to (20), thus 
getting more directly 

{Fip)-H{t))p=eP>*-^,-^^+b,+arp+a,y+--^H{t)= 

^ePit{bi+b.,t)H(t). 

The simplest way of arriving at this result is, of course, to 
apply the rule nr 3 of the table to the principal part of (20) 

and to interpret p as in the other terms. 


Our investigations have brought us asymptotic or exact expres- 
sions of the different constituents of / (<). As we see, every sing- 
Flp) 

ular point of — leads to an exponential factor, and obviously 


the term corresponding to the singularity to the extreme right 
will dominate over the other terms more and more as t tends 
to infinity, thus representing the ultimate form of the function 
/(/). If there are several such singularities, thus situated upon 
a line parallel to the imaginary axis, of course the sum of their 
contributions determines the final character of f(t). 


As an application of the general method developed we shall 
now derive the asymptotic expansion of the function 


f{t)=F{p).H{t) = --=ly..-H{t) (22) 

Vp*+1 


(Cf. pages 89 — 90). According to the uniformity-rule we have 


V p +w 


Here we have two branch-points, qi=i and q-j— — i, but no pole. 
To get the contribution from qi we first transform (22) by the 
shifting-rule, thus obtaining (Cf. (19)): 
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or 


Vp+2^ 

1 / \— 3 

/(<)=^_.e ^v'p•(l-i|) -Hit). 


(23) 


Then we expand about the origin according to the binomial 
formula : 




(24) 


2-4 -.(4A;) ® ^ ZJ^ 2 i - {4k-2) ^ ‘ 


■{4k) 

&= 1 Ar= 1 

Hence the asymptotic form of the term Iq^ looks: 

r 1 d 
“ 71' • 


■{ ,.“+! (25) 

A:=l 


A;- I 

Further we have: 


2} ’ 2-4 ■■■ (4A:— 2) ^ j ^ 






and generally, if n > 1, 


1 1 

2 sjnt t 


pn+L g (f ) =. ( - 1 )n ■ 111: i) i (26) 

2» V'TTf p 


Accordingly we get the following explicit expression for the func- 
tion Iq^ {t) (The factor H (/), of course, may be dropped) : 



m 


2 M Tit V^L> 2-4---(4A:) 

fc=l 


+.--2(-i> 

k=\ 


^ [l-3- 5--(4fe-3)? 1 _l_j 

2-4 - (4&— 2) 2*^-^ t^-^\ 


(27) 


For brevity we will write 


^^•M(<)+iB(<)] (28) 

Correspondently we obtain the contribution from the other 
branch-point g '2 = — i; obviously this must be the conjugate form 
of /gj, i. e. 

(28a) 


At last, combining the expressions Iq^ and Iq^, we arrive at 
the required asymptotic form of the Bessel function Jo{t), viz: 

(29) 


where, according to (27), the coefficients A{t) and J? (^) may be 
written : 




Jl-3-5--(4A;— l)p 1 


2k 


(8<)^’ 




B 




. [1-3-5- 1 




.... (30) 


If only the first term of A is taken into account, the develop- 
ment of course degenerates into the formula (252), derived by 
a more heuristical reasoning. 

We will not continue further the analysis of formula (29) but 
refer the reader to textbooks of Bessel functions and e. g. Funk- 
tionentafeln von Jahnke-Emde as to its managing in numerical 
calculus. 
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6. On the commutativity of the shifting integrator e““P. 

It follows from the definition integral of an integrational 
product that 

{p) H {t)==F {t-a). .. (1) 

As usual / (^) is assumed to be of class nr 1 at least. Especi- 
ally is 

Accordingly, the primary function of e'~^P has a point of dis- 
continuity at / = a; if a < 0 it does not vanish identically when 
^ < 0. Obviously this function is not of class 1, and there is no 
corresponding integral operator. 

Now we shall prove the following important proposition, applic- 
able in cases where occurs combined with integral operators: 
If a > 0, the commutative law 


^FX(p(t)-=Fx e~^P • cp (t) 
holds for any function cp{t) of class 1. 

The proof runs: 


I 

R. M. = j^/(0)+ J'f {t-r)dr^X(p{t-a) 


-f (0)-9P (< 


t 


{t—T) (p {t—a) dT. 


L. M. 


M.=^^-{Fxcp (<))=[^/ ( 0)-9 (0+ J' f' {t-v) cp 


t—a 


=f (Q)-(p {t — a)-\- J' i' {t—a — r') • 90 (t') • dr' 

0 

or, if we put T'==T“-~a, 

t 

L, M.~f (O)-^) (^ — a)+ J" f {t—T) (p{T — a) dT, 
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But now, a being positive, we have 9) (r — a) he 0 in the interval 
0 < T < a, and accordingly we may write as well 


L. M.==f f' {t — a)‘dr=li. M. 


Remark. In Chapter 6 we showed a certain parallelism between 
the symbolic method and the integrational calculus. Thus the 
exponential factor was seen to correspond to the unit-function 
H (t). The shifting integrator has its analogon in the ex- 
ponential factor in fact, F denoting the complex vector 

representing e. g. the harmonic function 

/ (l)=^/m 

we get f{t — a) represented by the vector For that reason 

the multiplying by may be characterized as an operation of 
phase-shifting. 


1 . Extended theory of integral operators. 

In Chapters 3 and 4 have been defined the four fundamental 
algebraical operations with ^ integral operators. Accordingly, we 
are in a position of forming arbitrary rational functions of integral 
operators, e. g. such a function As 

a F^+b HxK 
C’0^xW+d 

where the capital letters may denote operators, the others ordi- 
nary numbers. However, the procedure of division, being essenti- 
ally equivalent to the solution of a Volterra equation, will as a 
rule lead to not closed expressions. 

There may be reason for developing to some more extent the 
parallelism between the theory of basic integrators and that of 
integral operators. We wilt start with deducing the operators of 
the successive derivatives of a function. The operators of / (t), 
f (t) niay be denoted by jP, 'F, ** F. 
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Let (p (t) be any function and form the product 


y>(t)=FX9)(t)=/(0) 


I 

fp{t)+ ’’ 


(t — r)(p{T)dr. 


By differentiation: 


t 

j^/' (0) + f (t— f ) • dr j X 9) (<). 


But, by our manner of writing: 

t 

[/'(0)+y/"(<-Odrj = ’F, 


and accordingly we must have, (p{t) being arbitrary, 



The next step of differentialion gives: 

”■^=4 (4^-4 / (^)> follows 


dt \dt dt 


and it is easily proved by induction that 


(1) 




d" 

’ dr' 


d* 






dt 




which equality is formally equivalent to formula (20) of Chapter 2* 
From the identities: 

'FXcp(t) = ^[F-nO)]Xrp{t), 

"FXcp{t)=^[’F-fmx<p{t). 


we may conclude the following statement: 

From a function F X (p (t) to the function ^^'^Fxcf{t) one 
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will always come on diflerentiating a function with the initial 
value zero. 

This fact will be of importance in the applications because of 
the following commuiative law: 

The product 


1 

dt 


F=Fx 


d 

~dt’ 


i. e. is commutative, provided that it operates upon a function 
vanishing at the origin. 

For the proof we have to compare the two functions 
y’i{t)=^FXqi(t) and yh(t)—Fx^^(f){t). 


According to the foregoing: 

yji (t)=f (0) -9)' (t)+'FX(p{t). 

Further we get by integration by parts: 

t 

V^(<) = /(0)-9)'(t )+ — {r)dr 
0 

i 

=/ (0) • 9 (t) +[/'(<— r )- 90 (t)]'+ y " r (t — T) (p(T) dx 

. « 

=/ {0)-9p' {t) — f {t) (p (0)+’FX 9p (t). 

t 

Hence yh(t) — ifh(t)y only provided that f/)(0) = 0. 


Now we will apply the operator I = [^fdT] to both members 

0 

of (3). Because of the lemma just stated we get 

/X’F’X9>(0=Zx|[F-/(0)]X9>(0=|-/X[F-/(0)]X9)(t) 
=[F-/(0)]X9,(t). 

and, 9p(t) being arbitrary, this necessitates that 

IX’F=F-f(0). 


( 5 ) 
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By repeated use of this important rule we get 

7nx<">i?’=2?’- [/ (0)+/-/'(0)H (0)] (6) 

1 

which equality ought to be compared with (2). 


Formula (6) furnishes us with a means of getting by purely 
operationj^l calculus the solution of a linear differential equation 
with constant coefficients under given initial conditions. The 
procedure may be exemplified on the simple case: 


d^x 

di 






If X and F denote the integral operators of x {t) and / (^) 
respectively, the first equation may be written 

^'XxH {t)+2^^^XxH (t)+(o^-XxH {t)=^FxH {t). 

Obviously the law of uniqueness must subsist also for integral 
operators, and hence we get for X the “differentiaf* equation 

”X+2/?-’X+o>"-X-F (7) 

This may be solved by multiplying by P and then employing 
formula (6). Thus the equation will appear in the “algebraized” 
form: 

[l+2^-I+o)^-P]xX=PxF+u+(u+2^u)^I ( 8 ) 

The operational factor before X may be written: 

t 

L=l+2^I+o)^P=^+ y' (2/J+w=(<-T)}-drj (9) 

0 

The corresponding inverse operator L~‘ can always be obtained 
according to the general method (i. e. by solution of a Volterra 
equation), and we get on multiplying (8) by 

X=L--'xPxF+u-L-^+{v+2M-L-^XI- 
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Hence the explicit time function becomes 

X {t)=Xi (t)+X2 (/), 
xi {t)=^Lr^xPx f (f), 

X2 (t) = [u-\-{i)+2^u ) • /] X L~^ XH{t). 


( 10 ) 


This result should be compared with the integralional solution 
(83) in Chapter 6 of the same problem; the expressions agree 
formally, if p be interpreted as 


At last we are going to examine the inverse operation of 
FXcp{t) in the case of degeneration where /(O) — 0, thus 

t 


^^=[y/'(<-r)dr]. 


( 11 ) 


Let (p (f) be any function of class 1 and form 


v^(0=jP’x 


l 

'' it- 


T) (p[T) dr ( 12 ) 


Proposed that the primary function / {t) belongs to class 2 at 
least, we can differentiate both members of (12), thus getting 

t 

f(0) (p(t)+ J r {t—T) (p(r) dT=y>{t) 

0 * 

or 

'^X9>(0 = |v-{0 (13) 

Now, if f (0) ^ 0, we know the inverse operator of 'F to exist 
and may thus apply (’i^)“' to both sides of (13). Because i/) (0)=0, 
we obtain 

9(<)=(’F)-ix|v’(0=|c^’)-^xv(<) (14) 

Accordingly, when /(0)=0, we may define the inverse operator 
of F by the equality 


( 15 ) 
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This must always have the character of an integro-differential 
operator. 

If also f (0)=0, while f (0) ^ 0 and / (t) is of class 3 at least, 
we may take another step upon the same way, defining in (15) 
by the expression 

Hence it follows that 

and we may formulate the general rule: 

If the function /(t) belongs to class (n+1) at least and if 
f (0)=/' (0)= — =/<»—'> (0)=0, / <“>(0) ^ 0, the inverse operator of 
the integral operator F is defined by the formula 

F- = ((n)/r)-.X = (16) 

Accordingly, y) (t) being a sufficiently regular function, we are 
now in a position of solving the operational equation 

FXg){t)=-y)(t) 

even in cases where / {t) and some of its derivatives vanish at 
the origin. 

It may be of interest to identify the inverse operator of the 

fundamental integral operator I with the derivative operator 

For that purpose we have only to put /(0) = 0, f [t)=H (/), thus 
getting 


8. Note about the differentiation and the integration of a 
basic integrator with respect to a parameter. 

Let / {x, t) be a function of time and a parameter x (e. g. a 
space-coordinate) such that 


/ (», t), ^ and 


fm 


t)dx 
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have basic integrators. Then the following propositions are 
immediate consequences of the integral theorem: 


1. If 


t) 

~’e~P^dT, 


then 


£y/(a:,T)e-»<ir=y!’'<£ 

0 0 


( 1 ) 




2. If 


()X 

oo h 


J' jT J' j (iT, r) • {x, t) do ^ • e~^^dT, . . 


(2) 


then 


b b 

J F {x;p)‘H{i)\dx^^J' F{x;p) dx^H 


(t). 


It may be leaved to the reader to establish these relations. 

It is sufficient for the fulfilment of the conditions (1) and (2) 
that for the range of iT-values under consideration the respective 
time functions belong to class 1. Indeed, then the necessary 
uniform convergence of the infinite integrals is guaranteed, p 
being situated in a suitable half-plane. Accordingly, we may 
state the theorem: 

Any function / {x^ t) of class 1 may be differentiated or inte- 
grated with respect to the par^imeter x by differentiating resp. 
integrating the corresponding basic integrator F {x; p) before its 
application to H (t), provided that the ditferentiated resp. inte- 
grated time function belongs itself to class 1. 


9. Table of principal formulas.* 

'• 7/^)'" «•*=*-' 

* This table presents only such formulas as have been deduced in the text. 
It may be referred to the indices 4, 5 and especially 8 of the Bibliography for 
a more extensive collection of operational formulas. 


L 
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2 . 


p — a 








P^+flt 


“ .ff (<)=cos (oi-H (t). 


p"-+o)- 


H {t)—sin lot-H (t). 


(p+(i)- s\n (f+o>- cos (f « • / ^ I \ xr/^v 

P- ■H(t)^e-<»-sini(ot+cp)-H{t). 


l^P{0y2^Py V'ipr) 




[-■-m 


Hit),a>0 (218) 


p e “■''"’•ff (<)= - 

, 2 y'Tt t^t 


- — a • \Jp \ ■ — ®- 

i-c '' ■H{t)^-^-e 


p /V'p'^+w’* — p 


(tr\ (0 


■H{t)=Jx(ojt)-H(t), k>0... (233) 


V — a»\/p-+oj- 


(234, 241) 


■ J ry>{i~o )) • H {t—a)—Jo{(x) -y^P — a^)-H {t — a). 



^^7 (<)=/(<+«). 

O {p-\-a)-f <D (p) • e“*-/ {t). 

^ (p) • f • <I> (p+a) ■ e~°* 'f (t)- 
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